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ABSTRACT 


This  report  summarizes  the  work  performed  during  the  second  year  of  an 
experimental  and  theoretical  research  program  concerning  the  self-induced 
thermal  distortion  effects  on  cw  CO^  laser  radiation  propagating  in  absorbing 
media.  The  research  program  is  directed  toward  improving  the  understanding 
of  the  possible  limitations  imposed  by  the  atmosphere  on  the  propagation  of 
high  power  laser  radiation.  The  work  presented  is  concerned  with  essentially 
the  following  three  different  areas:  (l)  wind  simulation  experiments  using 
a  high  pressure  gas  cell;  (2)  studies  of  the  kinetic  cooling  of  a  gas  by 
absorbing  of  COg  laser  radiation;  and  (3)  studies  of  turbulence  effects  on  the 
self-induced  thermal  distortion  of  laser  beams. 
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SECTION  I 


INTRODUCTION 

The  United  Aircraft  Research  Laboratories  have  been  conducting  an  experimental 
and  theoretical  research  investigation  of  one  self-induced  thermal  distortion  effects 
on  cw  CO  laser  radiation  propagating  in  absorbing  media.  For  the  past  two  years 
this  research  has  been  under  the  joint  sponsorship  of  the  U.S.  Army  Electronics 
Command'  and  a  parallel  Corporate  sponsored  program.  The  research  is  directed  to¬ 
ward  improving  the  understanding  of  the  possible  limitations  associated  with  the 
propagation  of  high  power  C02  laser  radiation  in  the  atmosphere.  During  the  first 
year  of  this  program  considerable  progress  was  made  in  obtaining  an  understanding 
of  the  self -induced  thermal  distortion  effects  on  laser  beam  propagation  in  an 
absorbing  medium  in  the  presence  of  a  uniform  transverse  wind. (Ref.  1)  Experi¬ 
mental  and  theoretical  results  were  obtained  showing  the  thermal  distortion  to  he 
determined  by  the  nondime ns.lonal  distortion  parameter  N  and  the  Fresnel  number,, 
which  are  defined  in  terms  of  the  pathlength,  the  source  wavelength  and  beam  size, 
and  the  properties  of  the  propagation  medium.  Using  liquid  CSg  as  the  absorbing 
medium  experiments  simulating  a  wide  range  of  atmospheric  conditions  with  wind  were 
carried  out  for  values  of  the  distortion  parameter  N  up  to  ~100,  which  represents 
the  strong  nonlinear  interaction  regime  where  extremely  seyere  thermal  distortion 
and  intensity  reduction  effects  occur.  These  results  clearly  show  the  seriousness 
of  the  convective  thermal  lens  effects  and  through  the  use  of  the  parameter  N  the 
results  also  provide  a  means  for  obtaining  quantitative  estimates  for  the  atmospheric 
limitations  that  may  be  encountered  in  a  variety  of  ’nigh  power  cw  laser  applications. 

This  report  describes  the  work  carried  out  during  the  second  year  of  the 
program,  which  covers  the  period  from  30  April  I97I  through  29  April  1972.  During 
this  period  the  studies  of  thermally  self-induced  propagation  phenomena  have  pro¬ 
gressed  in  essentially  the  three  different  areas  of  (1)  high  pressure  gas  cell  wind 
simulation  experiments;  (2)  experimental  and  theoretical  studies  of  the  kinetic 
cooling  of  a  gas  by  absorption  of  C02  laser  radiation;  and  (3)  experimental  and 
theoretical  -tudies  of  the  effects  of  turbulence  on  the  self-induced  thermal  dis¬ 
tortion  of  iV-scr  radiation.  Section  II  describes  the  high  pressure  gas  cell 
experiments  ..a  which  results  obtained  for  the  thermal  lens  effects  of  wind  under 
large  N  conditions  are  examined  as  functions  of  the  focal  range  and  the  total  path 
attenuation.  In  addition  to  providing  data  for  estimating  the  atmospheric  limita¬ 
tions  on  high  power  propagation  for  a  wide  range  of  thermal  distortion  conditions, 
the  experimental  results  in  this  section  arc  also  important  for  providing  a  means  of 
checking  the  existing  nonlinear  propagation  codes  (Re'-.  2-6),  in  the  large  distortion 
regime  where  errors  and  instabilities  are  likely  to  become  important.  The  results 
of  the  focusing  studies  indicate  that  the  optimum  transmitter  focus  can  be  modified 
by  the  thermal  distortion;  and,  in  particular,  it  is  found  that  it  is  possible  for 
a  higher  target  intensity  to  be  achieved  wi«h  a  collimated  beam  than  with  a  focused 
beam  because  of  the  thermal  lens  effects.  This  suggests  the  technique  of  varying 
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the  transmitter  focal  range  as'  a  simple  first  step  toward  minimizing  the  thermal 
distortion  effects.  In  Appendix  C  &  related  experiment  is  described  in  which  the 
effects  of  astigmatism  on  the  source  beam  phase  distribution  are  examined  as  a 
possible  m^ans  for  reducing  the  thermal  distortion  by  a  wind. 

The  first  experimental  observation  of  the  kinewic  cooling  of  a  gas  has  been 
achieved  under  the  present  contract  using  an  electrically  pulsed  atmospheric  pressure 
COp  laser.  The  details  of  this  work  have  been  published  and  are  included  here  in 
Appendix  D.  In  Section  III  a  _.*ief  theoretical  treatment  of  the  kinetic  cooling 
effect  is  given  for  both  the  cases  of  pulsed  and  cw  laser  radiation.  From  these 
considerations.,  expressions  are  derived  for  the  magnitude  and  duration  of  the 
transient  density  variations.  Expressions  are  also  derived  for  the  steady-state 
cooling  effect  that  can  arise  with  a  cw  beam  in  a  transverse  wind;  and,  results 
obtained  using  them  in  the  nonlinear  propagation  code  are  included. 

Section  IV  describes  the  theoretical  and  experimental  efforts  to  determine, 
the  importance  of  atmospheric  turbulence  effects  for  the  propagation  of  high 
power  cw  CO^  laser  radiation.  A  theoretical  model  is  described  in  which  the 
effects  of  velocity  or  mechanical  turbulence  are  assumed  to  be  dominant  (in 
comparison  with  temperature  turbulence)  an"  u.re  characterized  by  an  eddy-diffusion 
coefficient.  Laboratory  experiments  with  irtifically  generated  velocity  turbu¬ 
lence  in  a  2  in  long  cell  are  also  described'.  The  results,  which  were  obtained 
for  the  conditions  of  strong  turbulent  diffusion  in  comparison  with  the  mean 
flow  velocity,  are  in  qualitative  agreement  with  the  theoretical  model.  The 
velocity  turbulence  provides  a  symmetric  mode  of  heat  transfer  which  tends  to 
reduce  the  thermal  distortion  effects  and  replace  the  asymmetric  bending  and 
spreading  by  the  mean  wind  velocity  with  a  symmetric  blooming.  The  importance 
of  the  turbulence  effects,  however,  is  proportional  to  the  gustiness,  i.e.,  the 
ratio  of  the  rms-to-the  mean  wind  velocity,  which  in  the  atmosphere  is  typically 
10-30$  or  less.  Thus,  mechanical  turbulence  effects  are  expected  to  represent 
at  most  only  a  small  perturbation  on  the  thermal  distortion  of  laser  beams  in  the 
atmosphere . 

Experimental  results  showing  a  comparison  of  the  distortion  with  beam  trans¬ 
lation  with  that  obtained  in  flowing  gas  in  the  50  cm  wind  tunnel  are  included  in 
Appendix  A.  These  results  indicate  that  the  poor  agreement  with  theory  of  some 
of  the  previous  wind  tunnel  data  (e.g.,  see  Ref.  1,  Fig.  27)  may  be  due  to  the 
thick  boundary  layers  associated  with  the  low  wind  velocities  rather  than  the 
effects  of  conduction  or  gas  heating  cited  earlier.  Thus,  because  of  the  possible 
interference  of  boundary  layer  and  turbulence  effects  it  is  clear  that  beam  trans¬ 
lation  experiments  are  preferred  over  the  use  of  small  wind  tunnels  for  simulating 
atmospheric  propagation  in  a  laminar  wind. 
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In  Appendix  B  the' 'correction  factors  to  account  for  focused  beam  propagation 
and  beam  sieving  in  the  thermal  distortion  parameter  N  are  derived.  The  previously 
derived  correction  factor  for  the  focused  beam  N(see  Ref.  1,  Eq.  (17))^  bas  been 
found  to  be  accurate  only  in  the  limit  of  small  amounts  of  focusing.  In  the  limit 
of  very  large  amounts  of  focusing  (i.e.,  when  the  ratio  of  the  initial  beam  size 
to  the  undistorted  beam  size  at  the  target  is  large)  the  correct  value  of  N  is  in¬ 
creased  by  approximately  a  factor  of  tvo  from  the  value  obtained  using  only  the 
initial-to-target  beam  size  ratio  (Ref.  1,  Eq.  (17))  to  account  for  the  focusing. 
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SECTION  II 

Large  N  Wind  Simulation  Experiments 

2.1  Introduction 

In  this  section,  experiments  using  a  high  pressure  gas  cell  are  described. 
Results  are  obtained  for  the  thermal  lens  effects  of  a  laminar  wind  under  condi¬ 
tions  where  the  distortion  parameter  N  becomes  large.  In  these  experiments  the 
thermal  distortion  is  examined  as  a  function  of  the  focal  range  of  the  source  ana 
the  total  path  attenuation.  The  advantage  of  using  a  high  pressure  gas  for 
obtaining  a  large  distortion  effect  with  a  wind  has  been  discussed  previously 
(Ref.  1);  and,  is  based  on  the  reduction  of  effects  due  to  thermal  conduction 
that  can  irterfeiewith  the  experiments  at  the  low  wind  velocities  which  are  used. 

In  addition,  difficulties  with  previous  gas  experiments  using  recirculating  wind 
tunnels  have  recently  been  shown  to  be  caused  by  the  boundary  layer  effects  as 
is  discussed  in  Appendix  A.  The  main  objectives  of  the  large  N  thermal  distortion 
experiments  with  the  high  pressure  gas  cell  were  to  examine  the  influence  of  vary¬ 
ing  at  and  also  to  study  more  carefully  rhe  effects  of  focusing.  Although 
large  N  thermal  distortion  results  have  been  obtained  using  liquid  CSg  ,  it  has 
not  been  possible  to  obtain  results  for  values  of  at  less  than  ~  2  due  to  the 
large  attenuation  of  CS^  at  10.6  p, . 

2.2  High  Pressure  Gas  Cell  Thermal  Distortion  Experiments 

The  experimental  arrangement  used  to  carry  out  the  thermal  distortion 
experiments  is  shown  in  Fig.  1.  The  experiments  consisted  essentially  of  propaga¬ 
ting  a  C0£  laser  beam  through  a  moving  cell  of  gas  oriented  vertically  (to  minimize 
natural  convection  effects)  which  simulates  the  effects  of  a  transverse  wind,  and 
then  monitoring  the  resulting  thermally  distorted  intensity  patterns  as  a  function 
of  the  laser  beam  power,.  The  laser  used  in  the  experiments  produced  approximately 
20  watts  in  the  TEM  mode  amd  was  stabilized  to  reduce  amplitude  fluctuation  in 
the  output  power.  ?n  the  first  series  of  experiments  to  be  described,  the  dis¬ 
tortion  effects  obtained  with  collimated  and  focused  beams  are  compared.  For  the 
focused  beam  case,  a  curved  mirror  (M^)  of  radius  2,15  m  was  used.  The  collimated 
beam  was  obtained  with  a  mirror  curvature  for  (M^)  of  6.78  m  radius.  The  high- 
pressure  gas  cell  consisted  of  a  5  cm  diameter  by  103  cm  long  tube  with  5mm  thick 
AR  coated  Ge  windows  at  both  ends.  For  all  the  experiments,  except  those  examing  it 

COg  at  a  pressure  of  150  psig  was  used  in  the 
nt  was  found  to  increase  from  ~1.8  x  10~^cnf  Vfc  t 
150  psig  or  ~  11  atm  pressure.  The  increase  in 
the  absorption  coefficient  is  due  to  the  effect  of  rotational  line  overlap  (Ref. 7) 
and  the  measured  increase  by  a  factor  of  ~  2.3  for  11  atm  is  in  I’easonable  agree¬ 
ment  with  the  value  of  2.6  that  has  been  reported  by  Christiansen,et .al.,  (Ref.  7). 
The  transmission  of  the  cell  windows  was  measured  to  be  ~  84  percent  with  Nr,at~5 
psig  in  the  cell.  The  absorption  coefficient  of  the  CO2  was  observed  to  increase 


at  dependence  .1  Section  2.4,  pure 
cell.  The  C0£  ibsorption  c^effjlcie 
atm  pressure  to  ~4.10  x  10  cm  at 
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somewhat  with  the  laser  beam  power  due  to  its  temperature  dependence  and  the 
influence  of  gas  heating.  For  example,  with  the  cell  stationary  the  CC>2  trans¬ 
mission  was  reduced  by  ~  7  percent  with  a  laser  beam  power  of  15  W.  With  the 
amount  of  cell  motion  involved  in  all  of  the  experiments  and  for  the  power  levels 
used,  however,  the  variation  of  the  CO2  transmission  with  power  could  not  be  observed 
within  the  accuracy  of  the  power  measurements  and  thus,  is  considered  to  be 
negligible . 

The  types  of  data  obtained  in  the  experiments  include  photographs  of  the 
distorted  laser  beam  patterns  obtained  using  Kalvar  film  and  intensity  profiles 
measured  by  moving  a  detector  equipped  with  a  35  micron  diameter  pin-hole  across 
the  beam  along  with  the  moving  gas  cell.  The  pertinent  parameters  for  the 
various  experiments  including  the  undistorted  laser  beam  radii  at  the  cell  input 
and  output  planes  together  with  the  cell  velocity,  absorption  coefficient,  path 
length,  etc.,  are  shown  in  Ihble  I  'with  the  experiments  identified  by  the 
appropriate  dates  and  focusing  conditions. 

In  Figs.  2,  3,  and  4  are  shown  the  data  obtained  with  the  collimated  beam. 

Fig.  2  contains  the  undistorted  intensity  patterns  and  profiles  along  the  wind 
direction,  i.e.,  the  direction  of  cell  motion.  This  data  is  shown  to  indicate 
the  deviation  from  the  ideal  symmetric  gaussian  beam  shape  which  is  assumed  in 
the  theory.  The  distorted  intensity  patterns  for  the  collimated  beam  are  shown 
in  Fig.  3  for  various  values  of  the  distortion  parameter  N,  where 

’  (-dn/dT)  2  P  z 
N  =  - 

.n  n  p  cpv  a~| 

Here,  n,  dn/dT,  p,  c  ,  v  and  a  are,  respectively,  the  refractive  index,  index 
change  with  respect  to  temperature,  density,  specific  heat  at  constant  pressure, 
velocity  and  absorption  coefficient  of  the  medium;  zis  the  absorbing  path  length, 
a.^  and  aQ  are,  respectively,  the  undistorted  laser  beam  radii  at  the  input  (source) 
and  output  (target)  planes,  and  P  is  the  power.  The  distortion  pai'ameter  N  is 
used  as  measure  for  the  strength  of  the  thermal  lens  effects  with  wind  and  was 
derived  originally  for  a  collimated  gaussian  beam.  As  shown  in  Appendix  B  the 
effect,  of  focusing  is  accounted  for  in  N  by  the  ratio  a^/a0  and  the  numerical 
correction  factor  q  (see  Fig.  B-l),  which  varies  from  1  to  2  as  ai/aQ  varies  from 
1  (collimated  beam  case)  to  ».  Although  this  form  of  the  focusing  correction  to 
N  applies  strictly  only  for  az  <  <  1,  it  is  also  used  for  finite  az  in  Eq.  (1) 
in  the  interests  of  keeping  the  expression  as  simple  as  possible.  The  values  for 
4  he  distortion  parameter  N  in  the  present  experiments  are  calculated  using  the 
measured  quantities  P,  aj_,  a0,  v,  at,  and  z  indicated  in  Table  I.  together  with 
the  known  properties  of  CO2  which  for  convenience  are  tabulated  in  Table  II. 

Here  t  i"  the  cell  length  and  to  account  for  the  small  air  path  (for  which  a  ~  0) 
between  bne  end  of  the  cell  and  the  detector  plane  we  distinguish  between  the  total 
path  z  and  t  aud  use  at  instead  of  az  in  the  term  inside  the  curly  brackets  in  Eq.(l). 
The  collimated  oeam  crescent  patterns  for  N=1.35,  2.7  and  5>^  in  Fig.  3  compare 
favorably  with  computer  results  obtained  by  Wallace  (see  Ref.  5,  Fig.  k)  for  the 
cases  z  =  a0  ,  1.5  z0  and  2  %0  ,  which  correspond  to  N  ~  1,  2.25  and  4,  since 
N  ~  (z/zQ)2  (Ref-  1»  P- 
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Quantitative  lata  is  obtained  for  the  thermal  distortion  in  terms  of  the 
magnitude  and  deflection  of  the  peak  intensity  as  found  from  the  measured  intensity 
profiles  shown  in  Fig.  4.  The  intensity  profiles  in  Fig.  4  are  grouped  together 
for  the  two  different  velocities  of  v  =  1  and  0.5  cm/s  to  show  that  the  inverse 
velocity  dependence  of  N  in  Eg*  '(X)  correctly  accounts  f or  the  observed  thermal 
distortion. 

In  Figs.  5-3,  data  obtained  for  the  focused  beam  case  are  shown.  The 
undistorted  beam-  oatterns  and  profiles  given  in  Fig.  5  indicate  that  the  beam 
deviates  somewnat  from  an  ideal  gaussian.  Comparison  of  the  distorted  intensity 
patterns  shown  Fig.  6  for  the  focused  beam  with  the  collimated  beam  data  in 
Fig.  3  shows  a  distinct  difference  in  the  nature  of  the  distortion  resulting 
from  the  convective  thermal  lens  effect.  In  particular,  the  small  islands  which 
form  in  the  distorted  patterns  for  the  focused  beam  at  the  larger  values  of  N  are 
absent  in  the  distorted  collimated  beam  patterns  which  retain  their  crescent 
shape  and  do  not  exhibit  the  tendency  to  break  up.  This  aspect  of  the  distortion 
has  also  been  found  in  numerical  calculations  of  focused  beam  propagation  as  is 
shown  in  Appendix  E.  In  Figs.  7  and  8,  the  proper  scaling  of  the  thermal  distor¬ 
tion  for  the  two  different  velocities  of  v  =  1  and  0.5  cm/s  is  again  clearly 
shown  by  the  measured  intensity  profiles. 

In  Fig.  9  a  series  of  distorted  patterns  and  intensity  profiles  are  shown 
which  indicate  progressively  the  influence  of  natural  convection  effects  as  the 
laser  beam  power  increases  from  P  =  7.2  W  to  14.4  W  for  the  wind  velocity  of 
0.5  cm/s.  The  effects  of  natural  convection  are  indicated  by  the  filling  in  of 
the  crescent  pattern.  An  estimate  of  the  natural  convection  velocity  for  a 
stationary  cell  oriented  vertically  can  be  obtained  from  the  expression  (Ref.  8) 


vnc~ 


gff  t  P 


2^t  cp  ix  T0 


(8) 


where  g  is  the  acceleration  of  gravity,  t  is  the  cell  length,  ^  is  the  viscosity  co¬ 
efficient  and  TQ  is  the3 ambient  temperature.  For  the  conditions  appropraite  to  Fig.  9, 
we  obtain  v  ~1.2xl0  where  P  is  in  watts.  Defining  the  characteristic  natural 
convection  time  Tnc=t/vnc  and  the  time  tw=  2  a^/v  associated  with  the  cell  volocity 

simulating  a  transverse  wind,  one  expects  the  ratio  t  /tw  to  be  near  unity  when 
natural  convection  effects  are  important.  Indeed,  for  the  power  levels  of 
~7  -  14  W  in  Fig.  9,  v  =  0.5  cm/s  and  a^  =  0.25  cm,  the  ratio  rnc/ tv  varies  from 
~3  to  2.  Evidently,  the  onset  of  natural  convection  is  associated  with  a  critical 
value  for  Tnc/rw  of  about  2.5  since  natural  convection  effects  appear  to  be 
unimportant  in  the  P  =  7.2  W  case  while  they  become  dominant  for  P  =  14.4  W  in 
Fig.  9. 


In  Figs.  10  and  11,  respectively,  the  normalized  peak  intensity  and  beam 
deflection  dependence  on  J5  are  shown  for  both  the  focused  and  collimated  beam 
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cases  c  The  normalized  peak  intensity  I__  is  the  ratio  of  the  distorted  peak 
intensity  to  that  intensity  which  would  De  obtained  in  the  absence  of  thermal 
distortion.  The  data  shown  in  Fig.  10  for  the  focused  beam  differ  considerably 
from  that  for  the  collimated  beam  foi  the  values  of  II  larger  than  three.  For 
the  focused  beam,  IpVr  decreases  with  II  at  a  rate  which  is  greater  than  l/ll,  for 

values  of  II  greater  than  II 3  -  The  rapid  decrease  of  Ij^r,  with  II  does  not 

occur  for  the  collimated  beam  until  II  becomes  greater  than  10.  Tne  difference 
between  the  collimated  and  focused  bean'  deflection  dependence  on  II  is  less 
pronounced  as  shown  in  Fig.  11.  However,  similar  to  the  normalized  peak  intensity, 
the  beam  deflection  for  a  collimated  beam  does  not  increase  as  rapidly  with  II  as 
it  does  in  the  focused  beam  case. 

A  different  way  of  presenting  data  for  the  effect  of  the  thermal  distortion 
on  the  laser  beam  intensity  is  shown  in  Fig.  IS.  -Here,  the  dependence  of  the 
peak  intensity  on  the  laser  beam  power  is  given.  At  very  small  power  levels,  the 
peak  intensity  increases  nearly  linearly  with  power  as  expected  for  propagation 
in  a  linear  medium.  However,  since  the  nonlinear  thermal  lens  strength  measured 
in  terms  of  N  also  increases  with  power,  a  maximum  in  peak  intensity  is  eventually 
reached.  Any  further  increase  in  the  laser  beam  power  Lends  to  reduce  the  peak 
intensity  because  of  the  severe  spreading  effects  produced  by  the  convective 
thermal  lens.  It  is  interesting  to  note  from  Fig.  12  that  the  maximum  target 
intensity  is  achieved  with  the  collimated  beam  rather  than  the  focused  beam,  but 
at  a  considerably  higher  value  of  laser  beam  power.  Also,  it  should  'be  noted 
that  doubling  the  velocity  from  0.5  to  1  cm/s  simply  doubles  the  power  level  at 
which  the  maximum  peak  intensity  occurs.  This  is  reasonable  since  the  maximum  of 
the  peak  intensity  versus  power  curve  should  always  occur  at  the  same  value  of 
II  =  Up  for  a  given  laser  beam. 

2.3  Liquid  CSg  Cell  Thermal  Distortion  Experiment 

For  convenience  in  canparing  the  high-pressure  gas  experiments  with  the 
liquid  CS2  experiments  and  also  because  the  latter  data  (Ref.  1)  has  been  replotted 
to  properly  account  for  the  effects  of  focusing  in  il  as  shown  in  Appendix  B,  the 
CS2  data  for  the  peak  intensity  and  beam  deflection  dependence  on  N  and  the  peak 
intensity  versus  power  are  shown  in  Figs.  13,  14,  and  15,  respectively. 

Comparison  of  these  results  with  Figs.  10  -  12  shows  that;  the  thermal  distortion 
results  follow  essentially  the  same  trends  in  both  the  liquid  and  nigh-pressure 
gas  cases.  It  is  also  found  from  Fig.  15  that  a  higher  peak  intensity  can  be 
achieved  at  a  target  with  a  collimated  beam  than  for  a  focused  beam  as  in  the 
high-pressure  gas  experiments.  As  in  Fig.  12,  the  label  i'lp  indicates  the  -value 
of  K  for  which  the  peak  target  intensity  reaches  a  maximum  with  increasing  power. 
Later,  it  will  be  shown  that  Up  corresponds  to  the  points  on  the  Ipg^  curves  in 
Figs.  10  and  13  where  the  slope  equals  minus  one. 
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2.4-  Dependence  of  Thermal  Distortion  onat 

With  the  high-pressure  gas  cell,  it  is  possible  to  vary  the  total  attenuation 
or,  equivalently  at,  by  varying  the  concentration  of  the  absorbing  gas  in  the  cell. 
To  examine  the  dependence  of  the  thermal  distortion  on  at  a  series  of  focused 
beam  experiments  with  three  different  absorbing  conditions  were  performed.  The 
focused  beam  was  obtained  using  the  Ge  lens  telescope  as  shown  in  Fig.  1.  The 
objective  in  these  experiments  was  to  determine  if  an  additional  dependence  on 
at  existed  beyond  that  implicit  in  the  expression  for  the  distortion  parameter  II 
given  in  Eq.  (l).  The  results  of  the  experiments  are  shown  in  Figs.  16-18;  The 
different  amounts  cf  total  cell  attenuation  were  obtained  by  using  a  mixture  of 
COp  and  nitrogen  to  give  at  =  0.13;  pure  C0o  for  which  at  =  0.44; and  a  mixture  of 
nitrogen  plus  a  small  amount  of  propylene  (C^IIg)  was  used  to  provide  at  =  2.3. 

In  all  cases,  the  total  pressure  in  the  cell  was  150  psig.  From  Figs.  17  and  18 
showing  the  dependence  of  the  normalised  peak  intensity  and  beam  deflection  on  N, 
it  is  clear  that  there  is  a  dependence  on  at  which  becomes  more  pronounced  with 
increasing  values  of  If.  Generally  speaking,  for  the  smallest  value  of  at,  the 
thermal  distortion  effects  are  more  pronounced  than  for  the  larger  value  of  at 
at  the  same  value  of  N.  This  is  borno-out  in  Fig.  16  where  the  distorted  intensity 
profiles  for  comparable  values  of  II  can  be  compared  for  the  three  values  of  at. 

Here  it  is  seen  tliat  for  at  =  0.13  the  thermal  distortion  is  more  severe  than  it 
is  for  at  =  0.44  and  2.3  at  the  same  value  of  II.  It  is  interesting  to  note  the 
axial  peaks  in  the  intensity  profiles  for  W  =  12  and  16  with  at  =  0.13,  which 
are  larger  than  the  peaks  associated  with  the  deflected  portion  of  the  beam. 
Although  ther^  appears  to  be  a  tendency  for  axial  peaks  to  be  formed  in  the 
intensity  profiles  fcr  the  larger  values  of  at  they  are  not  as  well  developed. 

A  possible  explanation  for  the  at  dependence  evident  in  Figs.  16-18  is  that  the 
values  of  N  given  by  Eq.  (1)  do  not  properly  account  for  the  effects  of  focusing 
when  az  is  large  (see  Appendix  B).  This  makes  sense  physically  since  it  is 
clear  that  focusing  will  have  little  effect  on  the  thermal  distortion  in  the 
limit;  large  az  because  the  intensity  is  reduced  by  attenuation  at  the  focal 
range  where  one  expects  the  focusing  to  have  the  greatest  effect  on  the  distortion. 

2.5  Focusing  Effects  on  Thermal  Distortion 

To  examine  the  effects  of  focusing  on  the  thermal  distortion  by  a  wind  the 
Ge  lens  telescope  was  adjusted  to  focus  the  beam  at  approximately  0.75  z,  z, 

1.25  z  and  at  some  range  between  1.25  z  and  «,  the  latter  being  the  collimated 
beam  case, where  z  is  the  range  to  the  target,  or  detector  plane.  The  differences 
between  the  thermal  distortion  obtained  with  those  focusing  conditions  are  shown 
qualitatively  in  Fig.  19 .  The  intensity  profiles  measured  as  a  function  of  laser 
beam  power  for  the  various  focusing  conditions  are  shown  in  Figs.  20  and  21. 

Values  of  II  were  not  computed  for  the  focusing  experiments  since  the  concept  of 
fl  is  of  questionable  value  for  beams  which  do  not  follow  the  usual  collimated 
or  focused  behavior.  For  this  case,  a  much  better  comparison  is  obtained  in 
t'-ms  of  the  peak  intensity  dependence  on  the  laser  beam  pow^r  shown  In  Fig,  22. 
Based  on  tne  data  shown  in  fig,  22,  the  maximum  target  intensity  is  achieved  with 
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the  bean  focused  at  the  target  or  when  F/z  =  1.  With  the  beam  focused  in  front 
of  the  detector  (f/z  =  0.7^-)  or  with  the  beam  focused  beyond  the  detector 
(F/z  =  1.26)  the  maximum  intensities  are  less  than  for  the  focused  case.  If  the 
beam  is  defocused  even  further,  i.e.,  F/z  >  1.26,  the  maximum  peak  intensity 
appears  to  be  slightly  higher;  and  presumably  for  the  collimated  beam  case,  the 
maximum  peak  intensity  will  be  greater  than  for  the  focused  beam  case  as  found 
earlier  in  Figs.  12  and  15 .  From  Fig.  22,  it  is  clear  that  the  thermal  distortion 
by  a  wind  changes  the  optimum  focal  range  for  the  laser  transmitter.  For 
example,  under  conditions  in  Fig.  22  where  the  laser  beam  power  is  4  watts,  a 
greater  target  intensity  can  be  achieved  by  focusing  25$  beyond  the  target  than 
by  focusing  at  the  target  range t  This  suggests  that  a  variable  transmitter  focus 
can  be  used  as  a  simple  means  of  reducing  the  thermal  distortion  effects  since 
under  extremely  large  distortion  conditions  (i.e.,  large  power  levels)  Fig.  22 
indicates  an  improvement  in  target  intensity  can  be  achieved  by  defocusing  the 
beam.  This  concept  is  related  to  experiments  discussed  in  Appendix  C,  in  which 
the  effects  of  astigmatism  are  examined  as  a  possible  means  for  reducing  the 
thermal  distortion  effects  of  wind. 


2.6  Summary  of  Wind  Simulation  Experiments 

The  simulation  of  the  tnermal  distortion  effects  of  a  laminar  wind  with  the 
liquid  CS2  and  high  pressure  gas  ceil  experiments  provides  results  which  can  be 
used  to  check  nonlinear  propagation  codes  and  to  develope  models  for  predicting 
the  limitations  imposed  by  thermal  blooming  on  various  high  power  applications. 

The  key  to  the  use  of  the  laboratory  simulation  results  lies  in  the  use  of  the 
dimensionless  parameters  N,  k  /z  and  az,  which  refer,  respectively,  to  the 
thermal  distortion  (i.e.,  the  strength  of  the  nonlinear  interaction);  the  Fresnel 
number  based  on  the  range  z  and  the  source  beam  radius  a^  ,  and  wave  number  k  - 
2rr/x  where  X  is  the  wavelength;  and  the  fractional  attenuation  factor  <yz.  By 
evaluating  these  parameters  for  the  situation  of  interest  one  can  then  refer  to 
the  normalized  peak  intensity  and  beam  deflection  results  in  Figs.  10, 11, 13; 14  or 
17,18  depending  on  which  are  most  appropriate  in  terms  of  the  degree  of  focusing 
and  the  total  path  attenuation. 

The  situation  can  best  be  s.jnmarized  as  follows.  The  distortion  parameter 
N  must  be  evaluated  using  Eq.  (1)  with  the  appropriate  values  for  the  properties 
of  the  medium  and  the  beam  radius  a.^  of  the  source,  which,  strictly  speaking,  is 
assumed  to  be  gauss ian,  and  the  beam  radius  a0  at  the  target  range  z.  Hie  degree 
of  . oc using  is,  of  course,  given  by  a^/a0  which  for  a  diffraction  limited  beam 
focused  at  the  range  z  is 'given  to  good  approximation  by  the  Fresnel  number 
k  a|  /z.  The  experimental  results  for  focused  beams  indicate  that  ,  which  is 

the  peak  intensity  normalized  by  the  undistorted  value  at  the  range  z,  drops  off 
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to  0.5  for  values  of  N  in  the  range  of  3  to  6  and  to  0.1  for  II  in  the  range  of 
9  to  18,  depending  on  the  magnitude  of  oz.  I™r  decreases  more  rapidly  "with  II 
when  ctz  «  1  than  when  yz  >  1  as  shown  in  Fig.  I?.  When  0*^  the  thermally 

self -induced  spreading,  tending  and  distortion  effects  are  probably  too  severe  to 
be  tolerated  for  most  applications  and  thus  N  should  probably  be  less  than  ~  10  -  20 
for  cases  of  practical  interest.  It  is  perhaps  worthwhile  mentioning  that  althoujn 
the  degree  of  focusing  has  been  limited  in  the  experiments  to  only  a  relatively 
small  range  the  results  mair  still  provide  a  reasonable  estimate  for  ...ore  strongly 
focused  beams  in  view  of  tne  uncertainties  associated  with  trie  effects  due  to 
turbulence  and  deviations  in  the  source  beam  profile  from  the  ideal  symmetric 
gaussian. 

For  the  case  of  a  collimated  or  diverging  beam  the  size  of  the  Fresnel  number 
is  important  since  it  determines  if  diffraction  effects  are  strong  enough  to  offset 
the  self-focusing  tendency  of  the  convective  thermal  lens  in  the  plane  parallel 
with  the  wind.  For  mosr,  situations  of  practical  interest  the  Fresnel  number  is 
expected  to  be  <  1C,  however,,  and  thus  the  self-induced  intensification  effect  will 

not  be  important  (Ref.  9) •  The  experimental  results  do  indicate,  however,  that 
j'REL  decreases  much  less  rapidly  with  N  for  a  collimated  beam  than  for  a  focused 
beam,  which,  no  doubt,  is  related  to  the  greater  spreading  and  the  formation  of 
the  central  lobe  in  the  distortion  of  the  focused  oeam.  It  is  important  to  note 
(see  Figs .  12,  15  and  22)  that  for  conditions  where  N  is  sufficiently  large,  it 
is  possible  for  the  target  intensity  to  be  greater  with  the  beam  collimated  (or 
perhaps  only  defocused  to  a  somewhat  greater  range)  than  when  it  is  focused  at 
the  Target.  This  can  occur  since  not  only  is  the  collimated  beam  N  smaller  by 
the  factor  a0/q  a^  than  the  value  for  the  focused  beam  but  also  IREL  decreases 
less  rapidly  with  N  for  the  collimated  beam  as  pointed  out  above. 


The  plots  of  IREL  versus  N  (see  Figs.  10,  13  and  17)  have  an  interesting  and 
useful  interpretation  if  one  considers  the  situation  where  only  the  laser  power  is 
variable  and  the  range,  wind  velocity,  absorption  coefficient  and  source  beam  size 
and  focus  remains  fixed.  In  this  case,  the  ordinate  on  the  1^^  versus  N  plots 
may  be  associated  with  the  power  since  II  is  proportional  to  the  laser  power  and 
the  remaining  conditions  are  assumed  constant.  The  peak  intensity  at  the  target 
range  z  is  given  by  I  (z)  =  Iu(z)  Irrl(N)  where  Iu(z)  is  the  undistorted  peak 

intensity  which  is  also  proportional  to  the  laser  power.  Clearly,  the  target  peak 
intensity  Ip  will  increase  with  power  so  long  as  IpjgL  does  not  decrease  with  N 
(or,  equivalently  p)  as  l/ll  or  faster.  From  the  Inversus  II  data  shown  in 
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Tigs.  10,  13  and  17  one  can  then  see  that  I  increases  with  the  laser  team  power 
until  the  slope  of  the  curve  equals  -1,  or  equivalently,  l/P*  Beyond  this 

point  where  3^  decreases  more  rapidly  than  l/u  the  value  of  Ip  actually  decreases 
with  increasing  power  because  of  the  influence  of  the  nonlinear  thermal  distortion 
effects.  -Thus,  the  value  of  H  for  which  the  log  1^,^  versus  log  II  plot  has  a  slope 
of  -1  corresponds  to  the  point  where  the  peak  target  intensity  reaches  a  maximum 
with  power  and  is  labeled  Up  in  Figs.  12  and  15.  The  value  of  Up  can  be  used  to 
determine  the  maximum  power  that  can  be  used  effectively  for  a  given  source  beam 
configuration  and  set  of  propagation  conditions..  This  power  and  the  associated 
value  of  Ipjgjj  then  determine  the  maximum  peak  target  intensity  obtainal.1  e  ur.der 
these  conditions . 

In  Appendix  E  numerical  results  obtained  by  Bradley  and  Herrmann  (Ref.  2) 
for  the  thermal  distortion  of  a  focused  beam  are  compared  with  a  distorted  laser 
beam  pattern  measured  in  a  liquid  CS2  experiment.  The  distortion  parameter  U 
is  used  to  make  t^e  comparison  between  the  simulation  experiment  and  the  computer 
calculation  which  is  for  high  power  propagation  conditions  in  the  atmosphere. 

The  good  agreement  of  the  results  establishes  the  validity  of  the  laboratory 
simulation  experiments  and  also  provides  evidence  for  the  usefulness  of  H  for 
scaling  the  thermal  lens  effects  of  wind  to  account  for  a  wide  range  of  conditions . 
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SECTION  HI 


KINETIC  COOLING 


3.1  Introduction 

The  kinetic  cooling  of  air  by  the  absorption  of  COg  laser  radiation  was  first 
discussed  by  Gerry  (Ref.  10)  and  has  since  been  studied  by  Glass  (Ref.  10)  and 
Wood,  Camac  and  Gerry  (Ref.  11).  The  cooling  effect,  which  is  transient  in  nature, 
causes  an  increase  in  the  gas  density  that  tends  xo  focus  the  laser  beam  and  thus 
the  effect  is  generally  regarded  as  being  beneficial  to  high  power  propagation. 

The  kinetic  cooling  effect  can  be  important  for  both  pulsed  and  ew  propagation. 

In  the  cw  case,  however,  a  transverse  wind  or  beam  motion  is  required  for  the 
kinetic  cooling  to  produce  a  steady-state  focusing  effect.  The  first  experimental 
evidence  for  the  kinetic  cooling  effect  (Ref.  12)  was  obtained  under  the  present 
Army  Contract  and  the  details  of  the  work  are  presented  in  Appendix  B.  The  kineti  ■ 
cooling  effects  were  observed  on  a  transient  basis  in  mixtures  of  COp  and  using 
an  electrically  pulsed  atmospheric  pressure  COq  laser.  Due  to  the  effects  of 
saturation  and  the  limited  sensitivity  of  the  interferometer  the  cooling  effects 
could  only  be  observed  with  CO2  concentrations  greater  than  5  percent.  As  a  result, 
the  characteristic  cooling  times  observed  were  limited  to  several  hundred  micro¬ 
seconds  as  compared  with  times  on  the  order  of  2-50  milliseconds  predicted  for 
normal  air  with  approximately  0.03  percent  COg  (Ref.  11).  Future  experiments  are 
planned  to  observe  the  kinetic  cooling  effect  in  normal  atmospheric  air  using 
either  a  nigher  energy  pulsed  or  a  high  power  cw  CO9  lasei-.- 

In  this  section  the  theory  for  the  kinetic  cooling  effect  in  air  is  briefly 
summarized  and  expressions  for  the  transient  gas  temperature  changes  resulting 
from  the  cooling  are  derived  for  the  cases  of  a  short  pulse  and  a  cw  C02  laser 
source  in  stationary  air.  These  expressions  are  needed  to  determine  the  optical 
effects  and  the  requirements  for  laboratory  experiments  investigating  the  kinetic 
cooling  process .  The  steady-state  temperature  expression  for  the  cooling  effect 
with  a  cw  beam  in  a  wind  or  with  beam  motion  is  also  derived  and  the  result  is  used 
to  modify  the  nonlinear  propagation  code  to  account  for  kinetic  cooling.  Numerical 
results  illustrating  the  kinetic  cooling  effect  with  wind  are  shown  in  Section  3.  . 

3.2  Summary  of  Theory 

Tiie  kinetic  cooling  effect  is  a  result  of  the  vibrational  energy  exchange 
processes  that  are  involved  ir.  the  absorption  of  10.6  y,  wavelength  radiation  by 
the  C02  in  air.  A  detailed  theoretical  treatment  of  the  absorption  process  has 
Led  to  a  simple  approximate  model  for  the  kinetic  cooling  in  air  (Ref, 11).  Accord¬ 
ing  to  this  model,  the  laser  beam  energy  absorbed  by  C02  is  stored  in  the  combined 
vibrational  levels  COg  (001)  and  IL  (v=l),  for  which  the  effective  collisional 
relaxation  time  is  t.  The  relaxation  time  for  the  lower  level  (c 0  (100)  )  of  the 
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absorbing  transition  is  essentially  instantaneous  in  comparison  with  t.  Thus, 
equilibrium  of  the  COg  (100)  level,  which  is  depleted  by  absorption,  is  rapidly 
restored  by  eollisional  energy  transfer  from  translation.  In  this  way,  the  gas  is 
cooled  until  the  stored  vibrational  energy  heats  the  gas  after  the  effective  re¬ 
laxation  time  t-. 

An  expression  for  the  gas  temperature  at  constant  pressure  is  found  from  the 
equation  for  energy  conservation  neglecting  heat  conduction  (Kef.  13) 


')*  o-1 


where  the  enthalpy  h  =  cpT  +  Ey  ,  I  is  the  laser  beam  intensity  and  p ,  cp  ,  T,  o-^ 
and  v  are  respectively,  the  density,  specific  heac  at  constant  pressure,  temperature, 


total  absorption  coefficient  and  velocity  of  the  gas. 


is  the  change  in  the  vibra¬ 


tional  energy  of  the  combined  CO2  (.001)  and  N2  (v~t)  levels  which  results  from  the 
absorption  of  10.6  p.  radiation  by  the  CO  ,  and  is  governed  by  the  equation 


Ev)  =  2-' 


«c  I(t)  -  p 


where  ct  is  the  CO^  contribution  to  the  total  absorption  coefficient  q'.j.  =  o,c+  cy^  ,, 
and  accounts  for  the  absorption  due  to  water  vapor.  The  energy  absorbed  by  water 
vapor  does  not  appear  in  Eq.  (4-),  since  it  is  rapidly  transferred  to  translation  and 
thus  leads  only  to  gas  heating  (Kef.  11).  The  factor  2.44  in  Eq.  (4)  is  the  ratio 
of  the  vibrational  energy  of  the  COg  (001)  level  to  the  10.6  p  photon  energy  hy, 
•which  accounts  for  the  fact  what  each  absorbed  10.6  p  photon  with  energy  E^^- 


results  in  a  vibrational  quantum  of  energy  E 
the  definition  for  h  and  Eq.  (4)  in  Eq.  (3)  ^ 

Ml  y .  vt  =  at 

a*  P  c 

v  p 

which  must  be  solved  together  with 


that  is  stored  in  the  N2  ,.  Using 


2.44  Zl' 


?Ey  = 


2.44  » 


l(t)  - 


*  Although  this  is  true  in  air  where  the  COg  concentration  is  typically  ~O.Q33  percent 
the  effect  of  the  lower  level  relaxation  time  cannot  be  neglected  when  the  COo  con¬ 
centration  is  >  10  percent  as  in  the  experiments  in  Appendix  D.  d 
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to  obtain  the  temperature. 

The  transient  temperature  changes  for  pulsed  and  cv  laser  radiation  are  now 
considered  with  the  gas  stationary  ( i.e . ,  v  =  0  in  Eas.  (5)  and  (6)  ).  For  the 
pulsed  case  the  pulse  width  T  is  assumed  to  be  much  shorter  than  the  relaxation  time 
t.  Thus,  for  l(t)  =  I_Tp6,(\t)f  where  6(t)  is  the  unit  impulse  function  and  I^t  is 
the  energy  density  of  the  pulse,  the  change  in  gas  temperature  following  the  laser 
pulse  is  given  by 


AT  (t) 


Oj.1 

top 


p  CT 


1 


2,44  — 


(7) 


If  or  /or  >  0.4l  transient  cooling  of  the  gas  will  occur.  The  cooling  effect  is 
greatest,  of  course,  if  (i.e.,  =0)  in  which  case  the  maximum  temperature 

decrease ,  which  occurs  at  t  =  o ,  is 


O'  I  T_ 

c  o  P 


|ATm|=  1.44  ~  ^  .  (Pulsed) 


(8) 


The  time  interval  for  which  the  gas  is  cooled  is  0.89  T  or  less  depending  on  whether 
o^/o?c  is  equal  to  or  greater  than  zero.  It  should  be  pointed  out  that  the  presence 
of  water  vapor  (i.e.,  or  ^  0)  will  not  only  reduce  the  magnitude  of  the  cooling  effect 
but  also  decrease  the  relaxation  time  t. 


For  the  case  of  a  cw  laser  beam,  we  assume  for  the  intensity  l(t.)  =  I  u(t^, 
where  u't)  is  the  unit  step  function,  and  from  Eqs.  ^5)  and  ( 6 )  the  temperature 
change  is 


AT^t)  = 


a .1 

t  o 


P  cT 


t 


2.44  t  (  l  -  e 

“t 


( 9 ) 


Assuming  =  aQ,  the  maximum  temperature  decrease  occurs  at  t  =  0.8-9  t,  and  its 
magnitude  is  given  by 


0.55 


O'  I  T 
CX0 

P 


(cw) 


(10') 
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Equations  ,(8)  and  (10)  can  be  used  to  determine  the:  requireaients  f  or  experimentally 
observing  t’ne  kinetic  cooling  effect  in  dry  air.  In  normal  atmospheric  air  the 
relaxation  time  t  is  reduced  by  the-  presence  of  0o  and  TV'  (Bef- s  11)  and  this  to¬ 
gether  "with  the  absorption  due  to  the  vater  vapor  vill  decrease  the  magnitude  of 
the  cooling  effect  as  indicated  by  Eqs .  (J)  and  (9) . 

For  the  case  of  a  ew  beam  vith  the  uniform  vind  velocity  v  if*  the-  x-direction 
a  steady-state  solution  for  the  gas  temperature  can  be  found  from  Eqs..  (5)  and  (6). 
Taking  the  Laplace  transform  of  Eq.  (6)  and  using  the  boundary  condition 
E^r  ;t  =  j)  =  0  leads  to  the  result 


,  2.44  an  fx 

E  (rjt)  = - 2_  / 

v  Vp  J 

*bc-vt 


-(x-x1  )/vt. 

I(x*,y,z)  e  dx*  .  (11) 


Substituting  this  expression  into  Eq.  (5)  and  setting  dT/ijt  =  0  while  letting 
t  -» 00  gives  the  equation  to  be  solved  for  the  steady-state  temperature 


dx  p  c  v 
P 


(1  -  2.44  _£) 


2.44  or  r 

‘  W-  / 

p  VV-T 


-(x-X*  )/vT  , 

I  (x',y,z)  e 


Integration  of  Eq.  (12)  from  to  x  gives  for  the  steady-state  temperature 


f(x)  -  T(_  «)  = 


7“  r  l(xiy^)  (1  -  6  e“(x-x  )/vr}  toi  ^  ^ 


vhsre 


6  =  2.44  —  - 


In  uhe  following  section  this  expression  is  used  to  modify  the  nonlinear  propagation 
code  to  account  for  the  kinetic  cooling  effect  for  a  cv  beam  in  a  vind. 
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3.3  Numerical  Results 

The  computer  code  used  to  describe  the  nonlinear  propagation  effects  for  a 
cw  beam  in  a  transverse  wind  is  described  in  Ref.  1,  Section  3.6.  To  account  for 
kinetic  cooling  in  the  numerical  procedure  the  appropriate  changes  must  be  made  in  the 
expressions  for  0  (z)  and  ^(z)  which  appear  in  Eqs., ''(36)  and  (39)>  respectively, 

of  Ref.  1.  Since’in  general 


8  (*) 

nt 


2  dn  7  m 
nQ  dT  t 


(15) 


where  V4.  is  the  transverse  gradient  operator,  the  effects  of  kinetic  cooling 
can  be  included  by  using  the  steady-state  temperature  expression  in  Eq.  (13)  to 
evaluate  Eq.  (.15)  with  the  result 


A  f  %  5 

x  (1-6)  I  +  _  f  I 

vr  ** 


-(x-x1  )/vt 
i  dx» 


A  r.X 

+  y  L  U 


Sy 


,  -(x-x')/vt.  a  , 

(1-6  e  '  "  )  dx\ 


where 


e  = 


(dn/dT)  01 1 

n  0  c  v 
oK  p 


(16) 


(17) 


Here, e  is  the  same  as  in  Ref.  1,  except  that  has  replaced  a.  Also, 

if  6  is  set  equal  to  zero,  or  equivalently  aQ  =  0  and  no  cooling  can  occur,  Eq.  (16) 
reduces  to  same  expression  as  in  Ref.  1,  Eq.'  (36),  as  required. 

The  expression  for  Y(z)  from  Ref.  1,  Eq.  (39)  is 

Y(z)  =  -s  (7t  +  !|1)  .  (a/2)  ,  (18) 

which  after  substitution  of  Eq.  (16)  and  performing  the  indicated  operations  be¬ 
comes 


L921004-8 


Y(z)  =  -1 —  \  2(1-6)  ^  +JL  I 

2  I  Bx  vt 


-  _L_  (1-  H  «  ,  P  I  e"(X-X')/TT  dx. 

(■vt)2  1  ““ 

+  r  aL  (i-6  .-(M,WW)  ax. 

+  i-  SI  f  a  (1-6  e-(x-x')/vT)  fc. 

rr  -  u  _ro  * 


(19) 


i  ay 


By 


Again,  if  6  is  set  equal  to  zero  in  Eq.  (19)  the  expression  reduces  to  the  proper 
result  for  pure  beaci  heating.  Equations  (16)  and  (19)  have  been  used  in  tne 
numerical  procedure  to  examine  the  convective  thermal  distortion  associated  with 
kinetic  cooling. 

The  two  parameters  6  =  2.44  ctc/o/^  and  (3  =  vt/&0  are  used  in  connection  with  the 
kinetic  cooling  to  characterize  the  nonlinear  propagation  effects  in  conjunction  with 
the  distortion  parameter  N  and  Fresnel  number  F  =  ka^  /z.  The  parameter  &,  which 
can  vary  from  zero  (pure  heating)  to  2.44  indicates  the  relative  importance  of  the 
transient  cooling  by  the  COg  absorption  as  compared  with  the  gas  heating  due  to 
the  absorption  by  water  vapor.  The  parameter  (3  is  the  ratio  of  the  relaxation  time 
t  to  the  gas  transit  time  a0/v,  where  aQ  is  the  l/e  intensity  radius  of  the  gaussian 
beam  source.  The  cooling  effect  will  be  most  important  when  t  is  comparable  to  or 
greater  than  a</v  (i.e.,  for  3  >  1),  which  means  that  the  vibrational  energy  due  to 

CC>2  absorption  that  is  stored  in  Ng  is  removed  from  the  vicinity  of  the  beam  before 
heating  the  gas  after  the  relaxation  time  t .  The  distortion  parameter  N  is  evaluated 
in  the  same  way  when  kinetic  cooling  is  included  as  for  the  case  of  pure  heating, 
with  the  total  absorption  coefficient  used  in  each  case.  With  kinetic  cooling  the 
parameters  6  and  3  must  be  specified  in  addition  to  the  parameters  N  and  F  which 
indicate,  respectively,  the  strength  of  the  nonlinear  interaction  and  the  import¬ 
ance  of  diffraction  effects . 

An  example  of  the  thermal  lens  effect  of  kinetic  cooling  with  a  wind  is  shown 
in  Fig.  23  for  a  collimated  gaussian  beam  with  6  =  2,44  (i.e.,  for  =  0)  and  3 
=  2.0.  The  constant  intensity  contours  representing  40  percent  of  the  peak 
are  shown  for  N  =  0.52  and  F  =  6.67  for  both  the  cooling  and  pure  heating  (i.e., 
6=0)  cases.  Comparison  with  the  40  percent  contour  for  the  undistorted  beam 
(i.e.,  for  N  =  0)  shows  that  a  strong  self-focusing  effect  is  associated  with 
the  cooling.  For  this  case  the  peak  intensity  is  increased  to  approximately 
three  times  the  undistorted  peak  intensity  and  the  beam  dimension  is  decreased 
roughly  by  a  factor  of  two.  It  should  also  be  noted  that  the  distorted  peak 
intensity  shifts  in  the  wind  direction  under  the  influence  of  cooling  rather  than 
into  the  wind,  as  is  the  case  when  beam  heating  occurs. 
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In  Fig,  2k,  the  relative  peak  intensities,  are  plotted  as  functions  of 

H  for  both  the  cooling  and  heating  cases.*  The  self -focusing  due  to  the  kinetic 
cooling,  vhich  for  the  case  plotted  is  for  the  conditions  §  =  2.44  and  3=2,  leads 
to  a.  more  rapid  change  in  the  distorted  peak  intensity  with  II  than  in  the  heating 
case,  where  both  focusing  and  blooming  effects  are  involved.  The  dashed  curves  in 
Fig.  2k  show  the  dependence  on  N  of  the  normalized  average  intensity 


_  2rra2 
■AVE  ~  Aeff 


=  2rra2 


Jl2  dA 

(F? 


(20) 


where  &eff  is  the  effective  beam  area  and  Jl  dA  is  the  total  power.  It  is 
interesting  to  note  that  although  IAVE  follows  IREr,  fairly  closely  for  the  pure 
heating  case,  there  is  a  greater  discrepancy  between  the  two  with  kinetic  cooling. 
The  reason  for  this  can  be  traced  to  the  beam  distortion  associated  with  cooling 
which  causes  the  central  portion  of  the  beam  to  contract  more  than  the  outer 
portions.  For  example,  with  the  N  =  0-52  case  in  Fig.  23,  the  diameter  of  the 
10  percent  contour  (not  shown)  is  reduced  by  less  than  a  factor  of  1.5  as  compared 
with  the  factor  of  two  decrease  in  the  diameter  of  the  40  percent  contour. 


The  normalized  beam  deflections  are  plotted  versus  il  in  Fig.  25  for 
both  the'  cooling  and  heating  cases.  For  small  values  of  N,  the  deflection  with 
cooling  is  ~3-</2,  which  is  equal  in  magnitude  to  the  heating  case  but  in  the 
opposite  direction.  As  N  becomes  larger,  the  deflection  with  cooling  increases 
less  rapidly  than  for  the  pure  heating  case. 

The  data  in  Figs  2^  and  25  represent  essentially  the  extremes  in  the 
behavior  of  the  magnitude  and  position  of  the  distorted  peak  intensity  with 
pure  heating  and  strong  cooling  conditions.  For  smaller  values  of  6  and  0  the 
deflection  and  Ipg^  curves  will  lie  somewhere  between  the  two  extremes.  In 
particular,  for  the  case  6  =  1.0,  where  the  absorption  due  to  CO2  is  ~  4l  percent 
of  the  total,  and  with  3  =2.0,  the  kinetic  cooling  effectively  balances  the 
heating  and  there  is  essentially  no  change  in  peak  intensity  or  deflection  as 
the  parameter  N  increases. 


*  The  quantities  F0  and  zQ  shown  in  Figs.  2h  and  25  refer,  respectively,  to  the 
Fresnel  number  and  range  for  which  N  =  1. 
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SECTION  I V 

TURBULENCE  EFFECTS  ON  THERMAL  BLOOf>iING 


4.1  Introduction 

Previous  studies  of  the  self-induced  thermal  distortion  of  high  power  cw 
CO2  laser  beams  have  ignored  the  possible  influence  of  atmospheric  turbulence. 

To  establish  a  realistic  model  for  assessing  nonlinear  propagation  effects  in 
the  atmosphere  it  is  necessary  to  be  able  to  estimate  how  and  to  what  extent  the 
turbulence  will  modify  the  thermal  lens  effects  from  those  predicted  on  the  basis 
of  a  laminar  wind.  The  temperature  fluctuations  associated  with  atmospheric 
turbulence  are  responsible  for  producing  the  so-called  linear  scattering  effects 
which  can  influence  both  low  and  high  power  propagation.  The  turbulent  velocity 
fluctuations,  on  the  other  hand,  have  no  direct  effect  on  low  power  propagation 
but  they  may  be  important  at  high  power  levels  because  of  their  interaction  with 
the  beam  heating  and  self-distortion  effects  due  to  the  absorption  by  COg  and 
water  vapor.  The  temperature  turbulence  effects  on  propagation  have  been 
studied  extensively  and  are  reasonably  well  understood.  Thus,  since  we  are 
mainly  interested  in  the  effect  of  the  turbulence  on  the  thermal  lens,  we  assume 
as  a  first  approximation  that  temperature  turbulence  effects  on  the  thermal  dis¬ 
tortion  may  be  accounted  for  by  an  appropriate  spreading  angle  for  the  mean  inten 
sity  profile.  The  effects  of  the  velocity  fluctuations  or  mechanical  turbulence 
on  the  thermal  distortion  thus  remain  and  are  the  main  sublet  of  concern  in  this 
section. 

The  principal  results  of  this  section  are  a  theoretical  model  accounting 
for  mechanical  turbulence  effects  and  laboratory  experiments  with  artifically 
generated  velocity  turbulence.  The  experimental  results  are  in  qualitative 
agreement  with  the  theory  and  show  that  velocity  turbulence  tends  to  reduce 
the  thermal  distortion  effects  and  replace  the  asymmetric  bending  and  spreading 
by  the  mean  wind  velocity  with  a  symmetric  blooming.  The  magnitude  of  the 
velocity  turbulence  effects  depends,  however,  on  the  ratio  of  the  rms  fluctua- 
tion-to-the  mean  velocity,  i.e.,  the  gustiness.  Since  in  the  atmosphere  the 
gustiness  is  typicaily  <,  10-30$,  it  appears  that  mechanical  turbulence  effects 
will  represent  at  most  only  a  small  perturbation  on  the  thermal  distortion  by 
the  mean  wind  velocity  in  the  atmosphere. 

In  Section  4.2,  the  theoretical  model  is  described  together  with  an 
approximate  solution  for  the  mean  intensity  profile  whir’i  is  used  to  show  how 
the  turbulence  modifies  the  thermal  self-distortion  effects  from  those  based  on 
a  laminar  wind.  The  laboratory  experiments  using  velocity  turbulence  generated 
in  a  2  long  cell  are  described  in  Section  4.3.  The  results  presented  in¬ 
clude  detailed  contour  measurements  of  the  distorted  COg  beam  patterns  obtained 
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in  cooperative  experiments  with  the  contract  monitor.  The  characteristics  of 
atmospheric  turbulence  and  its  expected  importance  for  high  power  propagation 
in  the  atmosphere  are  discussed  in  relation  to  the  laboratory  experiments  in 
Section  4.4. 

4.2  Theoretical  Model 

A  theoretical  model  to  account  for  tne  influence  of  velocity  turbulence  on 
the  thermal  self-distortion  of  a  laser  beam  by  a  wind  in  the  atmosphere  is  develop¬ 
ed  in  this  section.  Since  our  interest  is  mainly  in  the  effect  of  the  turbulence 
on  the  thermal  lens,  several  simplifying  assumptions  are  madei  First,  the  usual 
linear  scattering  effects  on  the  laser  beam  intensity  which  result  from  tempera¬ 
ture  turbulence  are  assumed  to  be  negligible  in  the  absence  of  absorption  or  beam 
heating.  This  is  a  reasonable  assumption  for  10. 6 n  wavelength  propagation  since 
these  effects  tend  to  be  small  for  paths  less  than  several  km.  Moreover,  under 
strong  scattering  conditions  due  to  normal  atmospheric  turbulence  the  eff  jts 
as  far  as  the  beam  heating  is  concerned  can  be  approximately  accounted  for  by 
assuming  an  appropriate  spreading  angle  for  the  mean  intensity  profile.  The 
second  assumption  is  that  the  absorption  Induced  beam  heating  is  sufficiently 
small  that  it  has  a  negligible  effect  on  the  structure  or  characteristics  associated 
with  the  normal  mechanical  turbulence  in  the  atmosphere.  Thus,  we  are  consider¬ 
ing  only  the  weak  interaction  limit,  which  is  reasonable  for  many  cases  of 
interest.  For  example,  with  a  100  kW  CO^  beam  1  M  in  diameter  propagating 
through  the  atmosphere  at  sea  level  with  a  1  mi/h  wind  the  absorption  due  to  C0? 
and  HgO  produces  a  temperature  rise  across  the  beam  of  order  &T  ~  0.02  °C,  which 
is  smaller  than  the  temperature  fluctuations  (typically  ~1  0 C)  normally  present 
in  the  atmosphere.  The  third  assumption  is  that  the  turbulence  is  homogeneous  and 
isotropic  and  can  be  characterised  by  a  scalar  eddy-diffusion  coefficient  K^. 

With  these  assumptions  for  the  theoretical  model  we  are,  in  effect,  saying  that 
any  turbulence  effects  on  the  thermal  blooming  are  predominatly  due  to  the  influence 
of  velocity  fluctuations,  i.e., mechanical  turbulence,  on  the  heat  transfer  process. 
The  problem  thus  reduces  to  finding  the  mean  temperature  profile  T  for  a  given 
turt”ience  structure  superimposed_on  a  mean  wind  velocity  transverse  to  the  laser 
beam  a”i«  The  mean  temperature  T  can,  in  principal,  then  be  used  to  calculate  a 
distorted  intensity  profile  using  the  same  techniques  as  with  a  laminar  wind.  It 
should  be  pointed  out,  however,  that  the  intensity  profile  obtained  in  this  way 
is  only  a  first  approximation  to  the  mean  intensity  which,  in  general,  will  also 
depend  on  the  higher  order  statistics  of  the  random  temperature  fluctuations. 

In  the  present  model  scattering  effects  which  can  arise  from  the  temperature 
fluctuations  due  to  the  interaction  of  the  turbulence  with  the  beam  heating  are 
also  ignored. 
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The  equation  for  the  mean  temperature  T  can  be  found  from  the  expression  for 
conservation  of  entropy  in  an  incompressible  fluid  by  the  usual  linearisation  pro¬ 
cedure.  Including  the  absorption  heat  source  term  a  f  (where  I  is  the  mean  intensity), 
ve  can  write  (Ref .  15) 


oT 

at 


U  ^ 

oy 


tT-va  =  X  t2  T 


O'! 

P~ c 
o  p 

-  A 


(21) 


where  T  =  T  -f  9  is  the  temperature*  and  v4x  +  uisthe  wind  velocity.  The  overbar 
indicates  an  average  and  9  and  u* are  the  random,  zero-mean  temperature  and  velocity 
fluctuations,  respectively,  (x  =  K/p0cD  is  the  thermal  diffusivity  of  the  gas )  .Equation 

21  differs  from  the  nonturbuLent  heat  equation  only  by  the  presence  of  the  term 

The  usual  procedure  at  this  point  is  to  make  use  of  the  fact  that  V-IT  =  0  to 
u-ve  =  rJ  •  "u$ ,  and  to  assume  that 


u-79. 

write 


-  -  (V  p 


a. 9 
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i i  ax, 


(22) 


..here  (K^)  is  a  second  rank  tensor  in  general.  The  coefficient  (Kh)i-  has  tne  same 
dimensions  as  the  Thermal  diffusivity  X  and  is  known  as  the  turbulence-or  eddy-trans- 
port  coefficient.  For  homogeneous  and  isotropic  turbulence  the  eddy-transport  co¬ 
efficient  becomes  a  scalar  constant  Kh  ,  and  for  the  steady-state  case  Eq.  (21)  now 
becomes 


u  . 
dx 


,  v  2-  o'! 

(K  +X)  V  T  =  _ 
h  Pc 

0  p 


(23) 


The  effect  of  turbulence  on  the  mean  temperature  profile  is  thus  accounted  for 
simply  by  a  diffusion-like  heat  transfer  term. 

An  order  of  magnitude  estimate  of  the  condition  necessary  for  the  turbulent 
diffusion  to  be  important  in  the  formation  of  the  thermal  lens  can  be  found  by 
comparing  the  heat  transfer  rates  l/T(3  and  1/tc,  associated,  respectively,  with 
the  diffusion  and  convection  terms  in  Eq.  (23).  Thus,  from  Eq.  (23)  we  obtain 
for  diffusion  1/t^  ~  +  Y)/a  and  for  the  convection  due  to  the  mean  wind 

velocity  1/tc  _  u/2a,  where  a  is  the  laser  beam  radius.  Since  for  air  the 
thermal  diffusivity  *  „  0.2  af/see  ,  which  can  usually  be  neglected  in  comparison 
with  the  eudy  diffusivity  Kh  ,  the  importance  of  the  turbulent  diffusion  in 
relation  to  the  mean  flow  is  proportional  to  the  ratio  T  /T  ~  4k  /Ua  The 

eddy  diffusion  coefficient  Kh  is  of  the  order  u't  (Ref?  jgf where  u* ‘is  the  rms 
velocity  fluctuation  and  *,q  is  a  characteristic  scale  length  for  the  turbulence. 
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I.ote  that  T  represents  the  deviation  from  the  neutral  vortical  temperature 
profile  vhicn  decreases  with  altitude  by  the  adiabatic  lapse  rate  of  l ()  0 
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If  the  diffusion  concept  is  to  be  applied  to  the  heat  transfer  from  the  laser 
beam  of  radius  a  by  the  velocity  turbulence,  the  eddy  diffusivity  can  only 
represent  turbulent  scale  sizes  X,  <  a.  Another  way  of  saying  it  is  that  the 
diffusion  process  is  a  valid  description  over  the  dimension  of  interest  only 
if  there  is  an  interaction  between  a-  large  number  of  eddies  ocurring  over  this 
dimension.  The  ratio  of  the  turbulent  diffusion  heat  transfer  to  that  by  the 
mean  flow  velocity  is  then  on  the  order  of 


(24) 


which  must  be  <  1  for  the  turbulence  to  be  important.  Since  in  the  atmosphere 
under  typical  conditions  u*/U  <  0.1  -  0.3,  we  do  not  expect  the  turbulence 
to  have  a  dominant  effect  over  the  thermal  distortion  associated  with  the  mean 
wind  velocity. 


According  to  Eq.  (23),  the  problem  of  determining  the  mean  temperature 
profile  for  the  laser  beam  in  an  absorbing  turbulent  flow  is  effectively  that  of 
combined  convection  and  conduction  heat  transfer.  Aleshkevich  and  Sukhorukov 
(Ref.  17)  have  considered  this  problem  and  obtained  the  solution  for  a  gauss iam 
beam  in  the  form  of  an  expansion  in  x  and  y  near  the  beam  axis.  Taking  the 
wind  to  be  in  the  x-direction  and  dropping  the  axial  heat  diffusion  term 
92T/  their  solution  for  'Eq.  (23)  can  be  written 


T(x,y, 


T(0,0) 


a  p 

o 

4  IT  K’ 


(25) 


where  K’  =  K^c^  +  K,  K  is  the  thermal  conductivity,  a  is  the  l/e  intensity 
radius  of  the  gahssian  laser  beam  and* 


T  = 

X 

ye2''"  [K(2va) 

1  V/ 

-  1^(2  y2)  ]  +  h  V"1  , 

(26) 

T 

XX 

=  [  (4  Y2  -  1  >  Kx 

2 

(2  y  )  -  4  y  K  (2  y  )  1  e"  Y  +  . 4 

(27) 

2  Y2 

T 

2 

2  y  2 

=  e  K  (2  y  ) 

1 

(28) 

yy 

2  Yd 

*  There  is  apparently  an  error  in  Ref.  17,  regarding  the  signs  of  T_,  and  T  . 
For  Eqs.  (4)  and  (5)  (in  the  same  reference)  to  be  consistent  with  Eq.  (2)  the 
minus  signs  in  Eq.  (2)  should  be  changed  to  plus  signs,  It  should  also  be  noted 
that  the  wind  is  taken  to  be  in  the  y-direction  in  Ref.  17 . 
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Here,  the  parameter  y  —  Us/4( K]j+  v )  is  the  ratio  oi  the  convecticn-to-eddy 
difxusion  heat  transfer  terms  in  Eq.  (23),  whirr  is  simply  the  reciprocal  of  the 
ratio  Tc/Td  defined  above.  The  quantities  and  T  ,  which  depend  only 

on  y,  are  given  in  Eos .  (26)  -•  (28)  in  terns  of  the  modified  Bessel  functions 
Kn*  (Ref  •  1^)  and  are  shown  in  Fig.  26  for  0  <  y  <  6.  Y7hen  y  -»  0,  diffusion 
heat  or  answer,  which  here  we  assume  to  De  due  to  the  turbulence,  becomes  dominant 
in  comparison  with  convection.  As  y  becomes  larger  than  one  the  convective  heat 
transfer  associated  with  the  mean  wind  velocity  becomes  dominant. 

It  should  be  pointed  out  that  the  expression  for  the  mean  temperature  in 
Eq.  (25)  is  only  accurate  for  |x/a|,|y/a|  «  1;  and,  in  addition,  Eq.  (25) 
is  not  suitable  for  use  in  the  iterative  computer  code  for  calculating  thermal 
distortion  effects  since  it  only  applies  for  a  gaussian  beam  while  the  computer 
procedure  requires  a  general  solution  for  arbitrary  intensity  profiles.  Equation 
25  can  be  used,  however,  to  derive  a  perturbation  type  of  expression  and  define 
a  scaling  parameter  to  account  for  turbulence  effects  on  the  thermal  distortion. 


Assuming  a  collimated  gaussian  beam  consistent  with  the  mean  temperature 
profile  in  Eq.  (25),  a  perturbation  expression  for  the  thermally  distorted  mean 
intensity  can  be  written  as  (assuming  &Z  «  l)  (Ref.  19) 


I(x,y,z)  =  I  (x,y,z)  e^Z  ^  f 


where  Iu  is  the  undistorted  gaussian  intensity  and 


<4*  =  - 


(dn/dT) 
2  n 


XJ  ^  'f  _  0  x  ^  _  o 

t  "  2  - 

a  Sx 


dy  ) 


(29) 


(30) 


which  accounts  for  the  thermal  lens  effects. 


Making  use  of  Eq.  (2?)  in  Eq.  (30)  the  expression  for  the  thermal  distortion 
exponent  $  is 

2 


=  -  IJY 


2  T  (Y)  /l)  +  T  (Y)  (1  -  2  *  ) 

x  \a/  xx  a2  ’ 


+  T  (Y)  (1  -  2  \ 


(3D 


*  In  Ref.  17,  the  are  referred  to  as  Macdonald  functions. 
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where 

(-dn/dT)  Q?  P  z2 
R  =  _  . 

n  n  P  c  U 
o  o  p 


(32) 


vhich  is  the  familiar  distortion  parameter  defined  for  the  nonturbulent  con¬ 
vective  thermal  lens  associated  vit'n  the  mean  velocity  U. 


It  can  be  seen  qualitatively  frc.n  Eq.  (31)  that  for  fixed  values  of  E  and 
mean  wind  velocity  U  the  effect  of  increasing  the  level  of  turbulence  as  measured 
by  will  lead  to  an  increase  in  on-axis  intensity  since  y  i/V 

intensity  profiles  through  beam  center  along  the  wind  direction  have  been 
calculated  using  Eq.  (29)  and  results  are  shown  in  Fig.  27  for  the  case  N  =  1. 

The  parameter  y-1  =  4iCu/Ua  ~  4u‘/U  provides  a  measure  of  the  relative  importance 
of  turbulent  heat  diffusion  and  since  in_Fig.  27  N  is  constant,  if  we  consider 
the  laser  beam  parameters  and  mean  wind  U  to  be  fi^ed,  the  changes  in  y  may  be 
regarded  as  changes  in  Kh  or  u* .  For  values  of  y  <1  there  is  little  effect 
on  the  Intensity  profiles  which  are  decreased  and  shifted  into  the  wind  character¬ 
istic  of  the  wind  or  convection  dominated  thermal  lens.  As  the  value  of  y~L  (or 
equivalently  Kh  or  u1 )  increases,  the  turbulent  diffusion  heat  transfer  reduces 
the  temperature  gradients  and  tends  to  restore  thg.  intensity  profile  to  the  un¬ 
distorted  condition  (N  =  0)  as  shown  in  Fig. (27). 

In  Fig.  23,  the  dependence  of  IR  ,  which  is  the  on-axis  mean  intensity 
normalized  b>  the  attenuated  undistorted  value,  on  the  parameter  y  ~  4u*/U  is 
shown  for  the  values  of  N  =  0.5  and  1.  As  the  turbulen£  diffusivity  Kh  (or 
equivalently  the  rms  velocity  u1)  increases,  so  that  y~  increases  from  ~  1  to 
10,  the  on-axis  intensity  is  approximately  doubled  for  the  case  N  =  1  and  for 
the  caxe  N  =  0.5  there  is  an  increase  by  a  factor  of  ~  1.4.  The  increase  in  on- 
axis  intensity  by  the  turbulent  diffusion  is  due  mainly  to  the  reduction  in 
beam  bending  which  can  be  seen  in  Fig.  27-  If  the  turbulent  diffusion  parameter 
y"1  can  continue  to  be  increased  beyond  10  the  thermal  lens  effects  can  eventually 
be  eliminated  as  sho^n  in. Fig.  28. 


*  It  should  be  pointed  out  that  since  the  temperature  expresssion  (Eq.  25) 
neglects  terms  higher  than  second  order  in  x  and  y,  the  intensity  profiles 
in  Fig.  27  are  only  accurate  for  |x/a|«  1;  and,  as  N  and  y  increase  the  curves 
become  increasingly  inaccurate  at  the  off-axis  points.  This  is  why  the  profiles 
do  not  converge  for  y-1  «  1  to  the  perturbation  results  shown  in  Eef .  19,  Fig.  2, 
as  they  should.  The  results  in  Fig.  27  are  correct  near  the  beam  axis,  however, 
and  are  also  useful  for  showing  qualitatively  the  effects  of  turbulence  on  the 
thermal  distortion. 
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4.3  Turbulence  Cell  Experiments 

The  approach  used  in  the  experimental  program  is  to  examine  the  thermal  self¬ 
distortion  of  a  laser  beam  in  a  laboratory  simulation,  using  artificially  generated 
velocity  turbulence.  One  of  the  considerations  involved  in  finding  a  suitable 
method  for  generating  the  turbulence  is  the  requirement  of  a  small  mean  flow  velocity 
to  enable  strong  thermal  lens  effects  to  be  obtained  with  the  small  30  W)  laser 
power  available.  Another  constraint  is  that  the  structure  of  the  turbulence  should 
include  a  range  of  scale  sizes  from  much  smaller  to  on  bhe  order  of  or  larger  than  the 
laser  beam  diameter  to  provide  a  significant  interaction  between  the  beam  and  the 
turbulence.  Although  the  use  of  grids  was  considered  initially  this  approach  was  not 
pursued  since  it  appeared  doubtful  that  fully  developed,  homogeneous  turbulence  would 
be  obtained  at  the  required  low  flow  velocities  and,  hence,  small  Reynolds  numbers  . 
After  investigating  several  types  of  nozzle  arrays  the  necessary  requirements  appeared 
to  be  met  in  a  turbulence  cell  consisting  of  a  series  of  opposing  jets  enclosed  in  a 
box  with  perforated  walls  on  either  side  as  shown  by  the  sketch  in  Fig.  29-  After 
initial  tests  with  a  3  in.  long  model  this  arrangement  was  extended  to  a  2  m  long 
cell  for  providing  a  region  of  nearly  homogeneous  and  isotropic  turbulence  with 
essentially  zero  mean  velocity.  The  effect  of  a  mean  wind  velocity  is  obtained  by 
moving  the  laser  beam  across  the  stationary  turbulence  cell  as  indicated  in  Fig.  29. 
Gelf-induced  thermal  distortion  experiments  are  conducted  using  pure  C02  in  the 
turbulence  cell  with  the  turbulence  level  being  monitored  by  measuring  the  gas 
pressure  at  one  of  the  plenum  chambers  feeding  the  tubes  with  a  water  manometer. 

The  turbulence  cell  initially  was  not  enclosed.  However,  it  was  found  necessary  to 
enclose  it  within  a  box  with  an  open  port  at  the  top  to  prevent  the  stratification 
and  mixing  of  CO2  with  air  which  produced  undesirable  refraction  or  bending  effects  at 
the  lower  turbulence  levels. 

In  Fig.  3Q  results  for  the  turbulence'  effects  on  a  stationary  beam  are  shown  with 
images  of  the  C02  beam  pattern  recorded  using  Kalvar  film.  With  no  turbulence  in  the 
cell,  which  was  oriented  horizontally,  a  vertical  natural  convection  velocity  is 
established  that  leads  to  the  downward  deflection  and  characteristic  kidney  or  crescent 
shape  seen  in  Fig.  30.  With  turbulence  in  the  cell  and  the  some  amount  of  absorption 
and  laser  power,  the  beam  patterns  become  symmetrical  and  decrease  in  area  as  the 
pressures  indicating  the  turbulence  levels  increase  from  1.5  to  16  cm  of  H2O.  This 
agrees  qualitatively  with  what  is  predicted  by  the  theoretical  model  in  Sec.  4.2. 

Figure  31  shows  photographs  of  the  patterns  near  the  exit  window  of  the  turbulence 
cell  of  a  HeWe  bcpi  that  is  superimposed  on  the  C02  beam.  The  HeNe  beam  was  well 
collimated  and  about  the  same  diameter  as  the  C02  beam  at  the  exit  window  of  the 
turbulence  cell.  The  HeNe  beam  patterns  in  Fig.  31  show  that  the  turbulence  effects 
on  the  stationary  CO2  beam  lead  to  the  symmetric  spreading  effects  that  are 
associated  with  conduction  or  diffusion  dominated  thermal  blooming.  The  nearly 
symmetrical  patterns  which  represent  ^  5  sec  exposures  indicate  that  the  turbulence 
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structure  is  stationary  and  nearly  homogeneous  and  isotropic  across  the  beam.  The 
pronounced  ring  structure  in  the  HeNe  beam  patterns  for  the  strong  thermal  blooming 
case  vith  27  W  in  the  COg  beam  are  believed  to  be  due  to  interference  effects  that 
are  enhanced  since  the  COg  beam  is  about  one-half  the  size  of  the  collimated  Helie 
beam  at  the  entrance  window  of  the  turbulence  cell.  Thus,  the  central  rays  of  the 
HeNs  beam  that  are  spread  by  blooming  eventually  intercept  and  interfere  vith  the 
outer  rays  that  are  unaffected  by  the,  COg  beam  induced  thermal  lens. 

Figure  32  shows  intensity  profiles  of  a  stationary  COg  beam  with  the  same  (6  cm  HgO) 
turbulence  level  for  different  laser  beam  powers .  A  35  p  pinhole  was  used  vith  the 
detector  which  was  scanned  slowly  (typically  1-3  minutes  in  duration)  across  the 
beam  to  provide  an  estimate  of  the  mean  intensity  profile  with  the  turbulence.  With 
ho  mean  wind  velocity  or  beam  motion  the  thermal  lens  effects  with  turbulence  should 
be  the  same  as  the  pure  conduction  blooming  case  but  with  the  thermal  diffusivity 
of  the  gas,  x  =  K/p0cp  replaced  by  the  turbulent  diffusion  coefficient,  K^.  For 
this  case  the  intensity  expression,  Eq.  (29),  becomes  (Ref.  20) 


I(r)  =  Iu(r)e"“z  expf-Nc[.2e~r2  - 1]} 


where  r  =  (x2  +  y2 


^u 


y 


(33) 


(34) 


and 


Nc  =  .(zdn/dT)  Pz 
TrnoPocp^ha 


1  -  P-JLaZl 

az 


(35) 


According  to  Eq.  (33)>  the  mean  intensity  profile  vith  a  stationary  beam  in 
turbulence  (or,  also  in  the  case  -yhere  thermal  conduction  is  dominant)  is  re¬ 
duced  near  the  axis  for  r  <  (tn2)2  ^  0.835; and,  for  r  >  0.835  the  intensity  is 
increased,  which,  as  the  thermal  blooming  becomes  stronger,  tends  to  change  the 
gaussian  profile  to  a  donut  shape.  The  changes  in  the  distorted  intensity 
profiles  vith  increasing  power  (which  for  fixed  turbulence  intensity  increases 
I’Jc)  are  seen  in  Fig.  32  to  be  in  qualitative  agreement  with  this  predicted  behavior . 
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Fie  mean  on-axis  intensity  normalized  by  the  attenuated  undistorted  value, 

1S  found  fram  Eg.  (33)  for  the  stationary  beam  in  turbulence  to  be 
given  by 


I„(o)  =  Iki_  =  e-"= 


I  e~az 


(36) 


Although  Eq.  (33)  is  a  perturbation  expression,  the  result  for  the  mean  on-axis 
intensity  in  Eq.  (36)  agrees  with  the  exact  analytic  expression  1^(0)  =  sech2N^ 

(Ref.  21)  within  ^  12$  for  E  <  1.  in  Fig.  33  measured  values  are  shown  for  the' 
mean  on-axis  intensity  as  a  function  of  power  at  the  turbulence  level  associated 
with  the  pressure  of  1.5  cm  of  HgO .  The  results  are  seen  to  follow  an  exporential 
decrease  as  predicted  on  the  basis  of  Eqs.  (35)  and  (36).  This  provides  a  con¬ 
venient  optical  method  for  determining  the  eddy- diffusion  coefficient  of  the 
turbulence  since  everything  in  the  expression  for  Kc  in  Eq.  (35)  is  known  except 
Kh  .  For  the  data  shown  in  Fig.  33  with  a  "turbulence  level"  of  1. 5  cm  of  H  0 
the  value  Kh  ~  3.5  cm2/sec  is  obtained.  This  value  can  be  compared  with  the2 
estimate  of  Kh  ~  1.2  cnr/sec  for  the  10  cm  H20  turbulence  level  which  is  based 
on  data  obtained  with  a  single  3/8  in.  long  hot  wire  probe.  The  rms  velocity 
fluctuation  measured  for  this  case  was  u*  ~  2  cm/sec  and  a  crude  estimate  for 
the  outer  scale  length  of  the  turbulence,  t+,  ~  0.6  cm  was  obtained  from  the  hot 
wire  anemometer  signal  with  the  probe  mechanically  scanned  at  25  cm/sec  through 
the  turbulence.  Additional  studies  comparing  values  of  K  obtained  optically 
with  hot  wire  anemometer  measurements  are  in  progress  to  check  the  theoretical 
model  and  establish  the  nature  of  the  turbulence  structure  as  a  function  of  the 
different  pressure  settings  given  in  cm  of  HgO.  The  pertinent  experimental 
parameters  for  the  stationary  beam  turbulence  results  of  Figs.  32  and  33  are 
given  in  Table  III. 

The  results  for  turbulence  effects  on  the  thermal  distortion  with  beam  motion  or, 
equivalently,  wind  are  shown  in  Figs.  34  and  35.  The  intensity  profiles  were  measured 
through  beam  center  along  the  direction  of  beam  motion  which  in  this  case  was  in  the 
downward  direction  (see  Fig.  29).  Several  traces  were  obtained  for  each  case  and 
the  results  indicate  both  the  repeatability  of  the  experiments  and  the  low  level  of 
intensity  fluctuations  produced  by  the  turbulence.  For  the  no  turbulence  cases  in 
the  lower  left-hand  corners  cf  Figs.  34  and  35,  the  peak  intensity  is  shown  to  be 
reduced  to  ~  36  percent  of  the  undistorted  value  and  shifted  into  the  wind  ~  1.65  a  . 
For  the  same  laser  power  and  beam  motion  conditions  but  with  velocity  turbulence, 
the  measured  intensity  profiles  on  the  right-hand  side  of  Figs .  34  and  35  show 
the  distortion  and  beam  deflection  being  reduced  with  the  on-axis  intensity  in¬ 
creased  as  the  pressure  setting  for  the  turbulence  increases  from  1.5  to  20  cm 
of  HpO.  The  only  important  difference  between  the  data  in  Figs.  34  and  35  is 
in  the  "wind"  or  beam  motion  velocities,  which  are  2  and  1  cm/sac,  respectively . 
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Since  the  power  level  used  in  Fig.  34  is  twice  that  in  Fig.  35 ,  the  same 
value  of  N  ~  11  is  obtained  for  each  case,  however,  as  is  also  indicated  by 
comparing  the  laminar  wind  (i.e.,  beam  motion  only)  cases  in  each  figure. 

For  the  maximum  turbulence  level  of  20  cm  H^O  with  U  2  cm/sec  (Fig.  3*0 
the  peak  intensity  is  increased  to  ~  60  percent  of  the  undistorted  value.  With 
U  =  1  cm/sec  (Fig.  35)  and  the  maximum  turbulence  level  of  1 6  cm  H2O  the  peak 
intensity  is  increased  to  ~  75  percent  of  the  undxstorted  value.  The  greater 
increase  in  the  peak  intensities  observed  at  all  turbulence  levels  with  U  =  1 
cm/sec  in  Fig.  35  as  compared  with  U  =  2  cm/sec  is  consistent  with  the  theoreti¬ 
cal  results  of  Figs.  2-7  and  28  since  y-1  -  4  K^/Us.  ~  **■  u'/u  for  35  is  twice 

that  for  Fig.  34  at  the  same  turbulence  level.  Eased  on  the  optically  determined 
value  of  Kh  ~  3.5  cm7seo  from  Fig.  33  for  the  1.5  cm  HpO  turbulence  level  and 
using_the  average  beam,  radius,  a  ~  0.26  cm  one  obtains  for  the  parameter  y~  = 

4  Kh/Ua  the  values  of  27  and  54,  respectively,  for  the  1,5  cm  Hr 0  case  in 
Figs.  34  and  35-  If  we  assume  that  K  ~  u'a,  the  rms  velocity  associated  with 
the  1.5  cm  of  H2O  turbulence  level  is  ~  13*5  cm/sec,  which  is  much  greater  than 
the  1  or  2  cm/sec  beam  motion  or  "wind"  velocities  used  in  Figs.  34  and  35- 


It  should  be  pointed  out  for  consideration  in  future  experiments  that  natural 
convection  effects  may  be  a  problem  under  experimental  conditions  similar  to  those  in 
Figs.  34  and  35  without  turbulence.  This  is  suggested  by  the  relatively  small  beam 
deflection  observed  for  the  laminar  wind  case  which  for  this  value  of  N  (~  11)  should 
be  about  %^4a0  instead  of  ~  1.65a0.  Indeed,  if  the  natural  convection  velocity  is 
estimated  (assuming  the  beam  to  be  stationary)  using  the  expression  (Ref.  20) 


vnc  ~ 


jot  g  P  \  l/3 


\PocpT , 


-3  1  2-2  -3 

with  the  conditions:  c  =  2  x  10  cm"  ,  g  =  9*8  x  10  cm  sec  ,  p  =  1.73  x  10 

g  cm""’,  Cp  =  0.86  J  g"1  °K'1,  T  =  290°K,  we  obtain  vnc  ~  1.66  P  1' ^  cm/sec,  with 

P  giver.  In  watts.  Thus,  for  P  =  8.3  V?  (Fig.  35)  vnc  ~  3.4  cm/sec  which,  since 

the  beam  motion  is  downward  and  at  a  velocity  of  only  1  cm/sec,  may  account  for 

the  smaller  beam  deflection  because  N  is  effectively  reduced  by  the  influence 

of  v  .  The  natural  convection  effects  may  be  reduced  by  orienting  the  cell 

vertically  or  possibly  by  making  the  beam  motion  in  the  horizontal  plane. 


More  detailed  results  for  the  thermal  distortion  effects  with  velocity 
turbulence  have  been  obtained  in  cooperative  experiments  with  the  contract  monitor. 
A  high  speed  (~  800C  rpm)  rotating  mirror  scanner  was  used  in  conjunction  with  a 
liquid  nitrogen  cooled  Hg  Cd  Te  detector  equipped  with  a  0.0039  in,  diam.  aperture 
to  obtain  time  resolved  (  ~  100  p,sec  scan  duration)  intensity  profiles  across 
the  beam.  The  beam  motion  across  the  scanner  was  provided  by  the  same  mirror 


28 


^S5B£sato(a 


L921004-8 


used  to  simulate  a  mean  wind  velocity  (see  Fig.  29).  The  intensity  profiles 
were  displayed  on  a  5^5  B  Tektronix  oscilloscope  with  a  1A7A  differential  pre¬ 
amplifier  and  were  recorded  on  35  M!-l  movie  film  at  36OQ  in/min.  Detailed 
constant  intensity  contours  were  obtained  from  the  measured  intensity  profiles 
approximately  39  cm  beyond  the  exit  window  of  the  2  ra  turbulence  cell  and  are 
shown  in  Figs.  3 6-46. 


The  undistorted  beam  contours  were  obtained  with  air  in  the  cell  (i.e ., 
no  absorption)  and  are  shown  in  Fig.  36.  The  input  and  output  1/e  beam  radii 
measured  from  individual  scans  through  beam  center,  perpendicular  to  the  direction 
of  beam  motion,  (the  beam  motion  is  from  right- to- left  in  Figs.  3 6-46)  are 
a^  =  0.21  cm  at.  the  cell  input  and  aQ  =  O.38  cm  in  the  output  or  detector  plane. 
The  contours  shown  in  Figs .  37-46  were  obtained  with  pure  CO^  in  the  cell  and 
based  on  a  transmission  measurement  with  F  ~  2  W  and  a  low  turbulence  level 
(  ~  6  cm  HgO)  the  value  <yt  =  0.37  iu  obtained  for  the  cell  absorption  coefficient- 
length  product.  The  laser  beam  power  entering  the  turbulence  cell  in  Figs.  37-40 
is  P  =  21  W  and  the  beam  translation  velocity  :s  /  =  2  cm/sec,  which  gives  the 
value  N  =  11  for  the  convection  dominated  distortion  parameter  (Eq.(l)).  The 
distorted  contours  are  shown  in  Fig.  37  with  the  COg  stationary  in  the  cell 
(i.e.,  the  laminar  wind  case).,  and  with  low  (6  cm  HgO),  medium  (10  cm  H^O)  and 
high  (16  cm  HgO)  turbulence  m  the  cell  in  Figs.  38,  39,  and  40,  respectively. 

In  Fig.  37,  with  no  turbulence  the  familiar  crescent  shaped  distortion  is  evident 
with  the  peak  intensity  (denoted  by  the  x)>  shifted  ~  2.6  aQ  in  the  direction  of 
beam  motion.  The  convective  thermal  distortion  spreads  the  beam  as  can  be  seen 
by  comparison  with  the  dashed  l/e  contour  obtained  with  air  in  the  cell,  and  the 
resulting  reduction  in  intensity  can  be  expressed  as 


I 

EFF 


AREA  l/e  CONTOUR-AIR 
AREA  l/e  CONTOUR. 


(38) 


For  the  no-turbulence  case  in  Fig.  37  Igpp  =  0.24.  With  velocity  turbulence 
superimposed  or.  the  same  2  cm/sec  beam  motion  the  distorted  beam  patterns  become 
more  symmetrical  with  the  peak  intensity  shifting  back  to  the  undistorted  beam 
axis  with  increasing  turbulence  in  Figs .  38-40.  Also,  in  Fig.  40,  with  high 
turbulence  1^  =  0.^75,  which  represents  approximately  a  55$  increase  as  a 
result  of  the  reduction  in  thermal  blooming  by  the  turbulent  diffusion. 


In  Figs.  4l  and  42,  the  cases  of  stationary  COg  and  medium  turbulence  are 

compared  for  the  conditions,  P  =  21  W  and  the  beam  translation  velocity  of  v  = 

5  cm/sec.  For  theoe  conditions  N  =  4.3,  and  for  the  no- turbulence  case  (Fig. 

the  peak  intensity  is  shifted  ~  2aQ  and  I^pp  =  0.33-  With  medium  turbulence 

(Fig.  42)  the  beam  deflection  is  reduced  and  the  crescent  shape  is  no  longer 

evident.  However,  no  increase  in  I__  is  obtained  with  the  turbulence. 

•  EFF 


1 

4i! 


In  Figs.  43-46,  the  cases  of  stationary  COg  and  low,  medium  and  high  tur¬ 
bulence  are  compared  for  the  conditions  of  P  =  11W  and  a  beam  translation  velocity 
of  5  cm/sec.  For  this  case  H  ~  2.2*5  and  with  no  turbulence  (Fig.  43)  the  peak 
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intensity  is  shifted  ~  1.75  aQ  and  the'  effective  intensity  is  =  0*54. 

With  low,  medium  or  high  turbulence  the  beam  deflection  is  reducea  as  expected  and 
the  thermal  distortion  tends  to  become  somewhat  more  symmetrical.  It  is  interest¬ 
ing  to  note,  however,  that  for  the  low  turbulence  case  (Fig.  44)  1^.,^,  is  reduced 


~  24$  below  the  stationary  CO^  case.  Although  I^_,  increases  to  about  the  same 
value  as  in  the  Lonturbulent  case  with  medium  turbulence  there  is  no  further 
increase  in  I  at  the  high  turbulence  level  (Fig.  46).  The  relative  ineffective¬ 
ness  of  the  turbulence  in  increasing  T  for  this  case  as  compared  with  the 
2  cm/sec  contours  in  Figs.  37-40  can  be  attributed  to  the  higher  beam  translation 
velocity  of  5  cm/sec,  which  for  the  same  turbulence  conditions  (i.e.,  the  same 
value  of  K^)  reduces  the  parameter  y"-*-  =  4  /Ua  by  a  factor  of  2.5. 


The  results  obtained  for  the  convective  thermal  distortion  in  the  presence 
of  velocity  turbulence  show  that  the  effect  of  the  turbulence  is  to  change  the 
asymmetric  wind  type  distortion,  which  is  a  combination  of  focussing  and  defocus- 
ing  distortion.  The  experimental  results  also  indicate,  in  qualitative  agreement 
with  the  theory  in  Section  h -2,  that  for  a  given  turbulence  level  its  influence 
on  the  thermal  distortion  is  inversely  proportional  to  the  mean  velocity;  and, 
that  the  rms  velocity  u1  must  be  much  larger  than  the  mean  U  for  any  appreciable 
reduction  in  the  thermal  blooming  to  result  from  the  turbulence  diffusion. 


4.4  Discussion 


Tiie  experimental  results  shown  in  Sec.  4.3  are  in  qualitative  agreement 
with  the  theoretical  model  developed  in  Sec.  4.2  in  which  an  eddy  diffusion 
coefficient  Kh  is  used  to  account  for  the  effects  of  mechanical  turbulence  on 
thermal  blooming.  Based  on  these  results,  it  is  found  that  velocity  turbulence 
in  the  presence  of  a  mean  wind  across  the  laser  beam  has  the  tendency  to  destroy 
the  combined  focusing  and  defocusing  of  the  pure  convective  thermal  lens  and 
replace  it  with  a  symmetric  two  dimensional  blooming  or  spreading.  Also,  accord¬ 
ing  to  the  theoretical  model,  the  velocity  turbulence  can  reduce  the  thermal 
distortion  and  spreading  effects  if  conditions  are  such  that  the  parameter 
y"1  =  4k, /Ua  is  much  larger  than  one.  In  the  laboratory  turbulence  experiments 
relatively  large  effects  on  the  thermal  distortion  were  observed  which  is  in 
agreement  with  the  preliminary  estimates  for  values  of  '»  1  that  are  based  on 
the  value  of  determined  optically  for  one  of  the  turb.’lence  levels.  Measure¬ 
ments  of  Kh  for  all  the  different  levels  of  cell  turbulence  used  in  the  experi¬ 
ments  are  la  progress  which  will  be  used  in  the  quantitative  evaluation  of  the 
thermal  distortion  results  in  terms  of  the  parameter  y~^. 

Although  the  cell  turbulence  structure  clearly  does  not  accurately  simulate 
that  found  in  the  atmosphere  the  laboratory  experiments  together  with  the  eddy 
diffusion  model  can  provide  a  qualitative  estimate  of  the  relative  importance 
and  effect  of  velocity  turbulence  on  atmospheric  propagation.  By  using  typical 
values  for  Kh  in  the  atmosphere  the  parameter  y-1  can  be  used  to  relate  the 
experimental  results  to  atmospheric  propagation  conditions.  The  eddy  diffusion 
coefficient  ~  u’t,  as  indicated  in  Sec.  4.2,  where  u*  is  referred  to  as  the 
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mixing  velocity  and  L  is  the  mixing  length.  Using  this  definition  for  K,  .  the 

*1  '**  Q  Jl  * 

parameter  y  ^  4(u'/u)  (-(,/a).  Since  the  length  t  can  be  >  lO^cm  corresponding 
to  the  largest  eddy  sizes  of  the  order  of  the  outer  scale  length  Lq,  and  4(u'/u) 
is  typically  of  order  unity,  it  -would  appear  that  can  indeed  be  »  1  for  the 
typical  beam  radius  a  ~  10  cm.  Thus,  on  this  basis,  our  original  predictions 
(Ref.  14)  -were  that  velocity  turbulence  in  the  atmosphere  could  be  expected  to  be 
an  important  factor  in  reducing  thermal  blooming  effects .  It  has  been  pointed 
out.  by  Bradley  and  Hermann  (Ref.  22),  however,  that  the  effects  of  eddy  sizes 
t  greater  than  the  beam  radius  a  cannot  be  represented  by  an  eddy  diffusion 
coefficient  since  the  diffusion  model  implies  the  interaction  of  many  eddies  over 
the  dimension  of  interest.  The  effect  of  the  larger  eddy  sizes  -l  >  a  must  be 
viewed  simply  as  a  slow  variation  in  the  wind  velocity  vector  that  can  be  accounted 
for  approximately  by  summing-the  squares  of  the  long  wavelength  components  with 
the  square  of  the  mean  velocity.  Thus,  if  we  assume  that  turbulent  eddy  sizes  -f, 
exist  that  range  from  much  smaller  to  much  larger  than  the  beam  radius  a,  which 
is  typical  for  atmospheric  propagation,  the  maximum  t  producing  a  diffusion 
effect  on  the  laser  beam  is  of  order  a,  the  beam  radius .  The  parameter  y-1  is 
thus  simply  of  the  order  4-  (u'/U),  which,  if  the  dependence  of  u'  on  l  is  ignored 

and  u'  is  taken  to  be  the  rms  velocity,  give**  y”1  <  1  for  typical  conditions  in 
the  atmosphere. 

A  somewhat  more  precise  estimate  of  y"1  for  typical  atmospheric  propagation 
conditions  may  perhaps  be  obtained  by  examining  published  data  for  K  .  Lettau 
(Ref.  23)^  for  example,  indicates  that  up  to  1  km  altitude  K  can  range  from  1 
to  Hr  cm  /sec  depending  on  the  mixing  length  l,  as  shown  in^Fig.  47.  Also  shown 
in  Fig.  4-7  is  the  molecular  diffusivity  x  ~  0.2  cm2/sec  for  air,  which  is  clearly 
unimportant  in  comparison  with  turbulent  diffusion.  If  we  assume  for  the  laser 
beam  radius  a  =  10  cm,  and  take  for  the  mixing  length  t.  =  a,  the  corresponding 
mixing  velocity  u'  varies  from  ~  1  to  10  cm/sec  as  shown  by  the  width  of  the 
shaded  region  in  Fig.  47.  Since  the  velocity  fluctuations  are  known  to  be  pro¬ 
portional  to  the  mean  velocity  0,  the  range  in  mean  velocity  associated  with  this 
variation  in  u'  is  given  approximately  by  the  mixing  velocities  of  the  largest 
eddy  sizes  which,  from  Fig.  47,  vary  from ~  10  to  200  cm/sec.  Using  these  extreme 
values  ^£or  U  and  taking  K  =  10  to  102  crn/sec  for  £  =  10  cm  in  Fig.  47,  we  find 
that  y  <  0.4  -  0.2  for  a  10  cm  radius  beam.  The  less  than  or  equal  condition 
accounts  for  the  fact  that  the  mean  velocity  must  either  be  larger  than  or  at 
least  as  great  as  the  largest  scale  length  velocity  fluctuations. 

If,  as  an  alternate  approach,  we  consider  the  Kolmogorov  similarity  theory 
of  turbulence  (Ref.  24)  we  find  that  except  for  the  smallest  eddy  sizes  of  order 
*‘Q  }  the  inner  scale  of  turbulence,  the  velocity  fluctuations  u'(^)  associated 
with  the  eddy  size  l  are  proportional  to  h1'  .  Taking  the  mean  velocity  U  as  an 
upper  limit  for  the  velocity  fluctuations  of  the  largest  eddy  sizes  Lq  ,we  then 
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have  that  u*  ~  U  (t/Lo)1^.  Since  l  must  be  <  a  for  the  eddy  diffusion  model 
to  be  valid  for  the  laser  beam  of  radius  a,  the  diffusion  coefficient  is  then 
=  u'(^)  i  <  u‘(a)  a,  and  thus  y"1  <  4  u'(a)/  U  ~  4  (a/L0)  >  vhich  for  the 

values  a  =■■  10  cm  and  L  _  =  lO^cm  gives  y~^<  4  (10)“^/^  ~  0.86.  Although 
this  gives  an  estimate  of  y  somewhat  larger  than  the  data  shown  in  Fig.  47,  the 
difference  is  not  significant  since  according  to  Figs.  ^7  ahd  28  the  turbulence 
effects'  on  the  thermal  distortion  are  negligible  for  y  <  1. 


From  these  considerations  it  appears  then  that  velocity  turbulence  effects 
on-  thermal  blooming  in  the  atmosphere  will  usually  be  insignificant.  This  con¬ 
clusion  is,  of  course,  in  agreement  with  the  intuitive  feeling  that  the  velocity 
fluctuations,  which  are  usually  smaller  than  the  mean,  cannot  be  as  effective  in 
removing  heat  from  the  beam  as  the  mean  flow  and,  hence,  have  proportionately 
less  effect  on  the  thermal  lens.  Possible  exceptions  when  velocity  turbulence 
may  become  important  are  when  the  mean  velocity  U  of  the  air  relative  to  the  beam 
approaches  zero  with  finite  turbulent  velocity  fluctuations  remaining.  Although 
it  is  doubtful  if  such  conditions  can  occur  with  a  stationary  beam  in  the 
atmosphere  except,  possibly,  under  strong  free  convection  mixing  conditions  near 
the  ground,  they  can  occur  with  beam  translation  or  slewing  when  the  beam  or 
points  along  the  beam  move  with  the  mean  flow  velocity.  In  this  case,  the  tur¬ 
bulent  diffusion  may  reduce  the  thermal  blooming  and  distortion  effects  from 
those  predicted  solely  on  the  basis  of  natural  convection  and  thermal  conduction 
in  the  regions  of  low  or  zero  relative  velocity. 
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TABLE  I 

EXPERIMENTAL  PARAMETERS 
HIGH  PRESSURE  GAS  CELL 

X  =  10.6  ji 

z  =  114  cm  (Total  Pathlength) 
t  =  103  cm  (Cell  Length) 

F  =  Focal  Range 


Bate 

?/z 

H'a  0 

a;  mrc 

v  cm/s 

ort 

6/23 

Focused 

1 

2.6 

3.1 

0.5,  1 

0.415 

7/6 

Collimated 

as 

0.9 

2.4 

0.5,  1 

0.415 

8/8 

Focused 

1 

7.4 

3.4 

1 

0.415 

8/8 

Collimated 

os 

o.97 

3.4 

1 

0.415 

S/9 

Focused 

0.74 

2.4 

3-6 

1 

0.415 

3/9 

Focused 

1 

7.5 

3-6 

1 

0.415 

8/17 

Focused 

0.74 

2.12 

3.4 

1 

0.415 

8/17 

Focused 

1 

7.1 

3.4 

1 

0.415 

8/17 

Focused 

1.26 

4.47 

3-4 

1 

0.415 

8/17 

Focused 

>1.26 

2.83 

3.4 

1 

0.415 

8/26 

Focused 

1 

7.1 

3-4 

1 

0.13,  0.44,  2.3 
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TABLE  ii 

MATERIAL  PROPERTIES 


cs2 

(293°K) 

Air  (1  atm) 
(300OK) 

COg  (1  atm) 
(300°K) 

n 

1.6319  (D) 

1.000263 

1.00417 

-dn/dT,  OK-1 

0.72  x  10“ 3 

0.875  x  10"6 

1.39  x  10"6 

p,  g  cm"3 

1.26 

1.175  x  10"3 

1.75  x  10" 3 

cp,  J/g°C 

0.95 

1.006 

0.86 

K,  W/cm°C 

1.6l  x  10-3 

2.6  x  10"^ 

1.66  x  10-1* 

Xj  cm2/sec 

1.35  x  10~3 

2.2  x  10"1 

1.08  x  10-1 

p,,  p,  poise 

3.66  x  103 

1.82  x  102 

1.46  x  102 
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TABLE  III 


EXPERIMENTAL  CONDITIONS  FOR  FIGURES  32  &  33 


STATIONARY  BEAM  IN  TURBULENCE 

Cell  Length  t  =  2  m 
Total  Bath  z  =  2.08  m 

a-t  =  0.371  (100$  C02) 
T  =  290°K 


Turbulence 

Level-CM  K20  mm  a0,  mm 


Fig.  32 

6 

1.23 

2.42 

Fig.  33 

1.5 

1.2 

1.9 
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FIG.  2 
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UNDISTORTED  FOCUSED  BEAM  DATA 
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THERMALLY  DISTORTED  INTENSITY  PROFILES 
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NATURAL  CONVECTION  EFFECTS 
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NORMALIZED  PEAK  INTENSITY  DEPENDENCE  ON  N 


HIGH  PRESSURE  GAS  CELL 
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NORMALIZED  BEAM  DEFLECTION  DEPENDENCE  ON  N 
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PEAK  INTENSITY  DEPENDENCE  ON  POWER 

HIGH  PRESSURE  GAS  CELL 
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NORMALIZED  PEAK  INTENSITY  DEPENDENCE 


C02  BEAM  IN  MOVING  CS2  CELL 
Q  *  —  2.2;  t  =  12*7  cm  v  =  1*5  cm/s«c 


N  X 


arbitrary  units 
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PEAK  INTENSITY  DEPENDENCE  ON  POWER 

COj  BEAM  IN  MOVING  CS2  CELL 
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FIG.  17 


INFLUENCE  OF  a  t  ON  THE 
NORMALIZED  PEAK  INTENSITY  DEPENDENCE  ON  N 


jTi] 


EFFECT  OF  FOCAL  RANGE  VARIATION  ON  INTENSITY  PROFILES - 

8  17  71 

HORIZONTAL  SCALE  =■  1.9  m  in  div 
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FIG.  20 


INTENSITY  DEPENDENCE  ON  POWER 
HIGH  PRESSURE  GAS  CELL 
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COMPARISON  OF  BEAM  DISTORTION  WITH  AND  WITHOUT  KINETIC  COOLING 
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EXPERIMENTAL  ARRANGEMENT 


L921004-8 


FIG 


SCANNING  DETECTOR 


TURBULENCE  EFFECTS  ON  STATIONARY  C02  BEAM 
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FIG.  30 
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HeNe  BEAM  SUPERIMPOSED  ON  C02  BEAM 
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FIG.  31 
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NORMALIZED  MEAN  ON-AXIS  INTENSITY  VS.  POWER 
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TURBULENCE  LEVEL  :  1.5  CM  HjO 
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FIG.  35 
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FIG.  37 
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MEASURED  CONSTANT  INTENSITY  CONTOURS 


BEAM  MOTION 
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FIG.  39 
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FIG.  40 
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FIG.  42 


MEASURED  CONSTANT  INTENSITY  CONTOURS 


BEAM  MOTION 
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FIG.  43 


MEASURED  CONSTANT  INTENSITY  CONTOURS 
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FIG.  44 
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FIG.  45 


MEASURED  CONSTANT  INTENSITY  CONTOURS 


BEAM  MOTION 
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FIG.  46 
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TURBULENT  DIFFUSION  DIAGRAM 

(FROM  H.  LETTAU,  “COMPENDIUM  OF  METEOROLOGY,”  1951) 
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APPEHDIX  A 

Comparison  of  Beam  Translation  with 
Wind  Tunnel  Flow 


As  a  result  of  some  previous  discussions  regarding  the  poor  agreement  wren 
theory  of  the  beam  deflection  results  obtained  in  gases  (see  Ref.  1,  Fig.  27), 
it  was  decided  to  compare  the  thermal  distortion  effects  obtained  with  beam 
translation  in  a  quiet  gas  with  those  obtained  using  the  gas  flow  in  the  cir¬ 
culating,  50  cm  path  wind  tunnel.  The  wind  tunnel  was  filled  with  a  mixture 
of  air  and  propylene  which  provided  enough  absorption  to  give  a  value  of  N  ~  k 
with  a  5  cm/sec  wind  or  translation  velocity  and  a  10  watt  laser  beam  power.  The 
distorted  profiles  were  measured  first  with  the  laser  beam  translating  across  the 
quiet  gas  and  then  with  the  laser  beam  stationary  in  the  wind  tunnel  gas  flow. 

The  hot  wire  anemometer  used  to  measure  the  wind  tunnel  flow  velocity  was  calibrated 
by  placing  the  probe  on  the  mechanical  scanner  used  for  beam  translation.  In  Fig.  A-l, 
the  results  are  shown  comparing  the  thermal  distorted  intensity  profiles  as  obtained 
by  the  two  different  methods.  There  is  a  significant  difference  between  the  distortion 
obtained  with  beam  translation  at  5  cm/sec  and  the  5  cm/sec  wind  tunnel  gas  flow 
velocity.  With  beam  translation,  the  shape  and  deflection  of  the  distorted 
intensity  profile  appear  to  be  in  reasonable  agreement  with  theoretical  predic¬ 
tion?.  In  the  wind  tunnel  flow  case,  however,  the  distorted  intensity  profile 
has  two  peaks  and  is  broader  than  with  beam  scanning  at  5  cm/sec.  Although 
increasing  the  wind  velocity  and  laser  beam  power  by  a  factor  of  two  leaves  the 
distortion  parameter,  H,  unchanged,  the  resulting  distorted  intensity  profile  in 
Fig.A-lis  clearly  seen  to  more  closely  approximate  the  distortion  obtained  using 
beam  translation.  Increasing  the  wind  tunnel  velocity  further  decreases  the 
the  distortion  parameter,  N,  while  improving  the  agreement  of  the  distorted 
intensity  profiles  with  theory.  Thus,  the  discrepancy  between  the  case  of  beam 
translation  and  the  use  of  wind  tunnel  gas  flow  appears  to  be  associated  with 
the  magnitude  of  the  flow  velocity. 

To  investigate  the  wall  boundary  layer  as  a  possible  cause  for  the 
observed  discrepancy,  the  wind  velocity  profiles  near  the  wall  were  measured 
for  the  two  free-stream  velocities  of  5  and  20  cm/sec.  The  results,  shown  in 
Fig. A-2,  indicate  an  increase  in  the  boundary  layer  thickness  from  *-  to  1;  inches 
as  the  velocity  decreases  from  20  to  5  cm/sec.  This  suggests  that  the  reduced 
velocity  associated  with  the  boundary  layer  near  the  wind  tunnel  entrance  window 
may  be  responsible  for  the  difference  between  the  gas  flow  and  beam  translation 
cases . 
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Turbulence  effects  do  not  appear  to  be  important  in  these  wind  tunnel 
gas  flow  experiments  since  the  distorted  profiles  are  observed  to  be  very 
repeatable  (e.g.,  the  data  in  Fig.A-1,  represent several  scans  of  the  distorted 
beam)  with  essentially  no  random  fluctuations . 

Based  on  the  results:  of  this  experiment,  the  relatively  poor  agreement 
wi th  theory  of  some  of  the  previous  data  obtained  using  wind  tunnel  flow 
(e.g.,  see  Ref.  1,  fig.  27)  may  be  attributed  to  the  boundary  layer  thickness 
rather  than  the  effects  of  conduction  or  gas  heating  as  cited  earlier.  Thus,  for 
the  largest  values  of  N,  where  the  largest  deviations  from  theory  were  observed, 
the  wind  tunnel  velocities  were  becoming  very  small  with  an  accompanying  increase 
in  the  boundary  layer  dimension,. 

In  summary,  it  is  now  clear  that  thermal  distortion  experiments  using 
wind  tunnel  gas  flow  must  be  carefully  monitored  to  avoid  boundary  layer  and 
possibly  turbulence  effects.  It  also  appears  that  with  beam  translation,  the 
experiments  can  be  carried  out  much  more  simply  and  with  accurately  simulated 
conditions  for  large  values  of  N,  even  in  gases. 
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APPENDIX  B 


The  Distortion  Parameter  for  Focused  and 
Diverging  Eeams  and  for  Beam  Slewing 


The  form  of  the  distortion  parameter  N  derived  previously  (Ref.  1,  ,Sec.  3.3) 
to  account  for  focused  beam  propagation'  has  been  found  to  be  correct  only  in  the 
limit  of  a  very  small  amount  of  focusing  or  divergence  of  the  beam.  The  general 
form  of  the  distortion  parameter  for  focused  or  diverging  gaussian  beams  in  a 
uniform  laminar  wind  is  derived  in  this  appendix.  The  form  of  N  appropriate 
for  the  slewing  of  a  collimated  gaussian  beam  is  also  given. 

The  change  in  intensity  of  a  laser  beam  due  to  gradients  in  the  refractive 
index  can  be  written  (n*-?.  19) 

=  exp  (Y)  (B-l) 

I  e-*2 
u 

where 

Y  =  "  S  (V  +  )  *  J  ZfL  dz->  dz*  ,  (B-2) 

o  t  I  on 

o 

Iu  is  the  indistorted  intensity,  i.e.,  the  intensity  in  a  homogeneous  loss-free 
medium,  n  =  nQ  +  §n  is  the  refractive  index  with  1 6 n/nQ |  assumed  to  be  «1, 

0?  is  the  linear  attenuation  coefficient,  z  is  the  propagation  direction  and  V 
is  the  transverse  gradient  operator.  For  the  case  of  a  convection  dominated  ^ 
thermal  lens  (i.e.,  a  transverse  wind  or  beam  motion  from  slewing  or  a  moving 
source)-  the  normalized  refractive  index  gradient  is  (Ref.  19)* 


^L. - =  e(z)  /x  I  +  y  x  dx'l 

n0  \  -  5y  | 

where 

(dn/dT)  a 

e  (z)  - -  , 

n0  p  cp  v 


(B-3) 


(B-4) 


^Although  kinetic  cooling  is  not  included  here  to  simplify  the  derivation  it  could 
easily  be  incorporated  in  Eq.  (B-3)  by  comparison  with  Eq.  (15)  and  (16).  The 
expression  for  N  is  the  same  in  both  cases,  however. 
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which  is  written  as  a  function  of  z  to  account  for  changes  in  the  properties  of 
the  medium  along  the  propagation,  path.  For  example,  the  altitude  dependence 
of  the  absorption  coefficient  and  the  properties  of  the  air  could  be  accounted 
for  along  a  non-horizontal  propagation  path;  or,  the  wind  velocity  may  be  non- 
uniform  along  the  propagation  path.  The  distortion  parameter  H  is  essentially 
the  dimensionless  factor  that  multiplies  the  geometric  terms  in  Eq.  (B-2)  if  a 
perturbation  expression  is  obtained  by  assuming  I  »  1  e  .  Here  we  wish  to 
obtain  the  general  form  for  ¥  that  applies  for  an  uncollimated  gauss ian  beam  with 
a  nonuniform  transverse  wind  or  beam  motion,  e^g.,  such  as  slewing  along  a 
horizontal  propagation  path. 


Thus*,  assuming  for  the  undistorted  beam  the  uncollimated  gauss ian  intensity 
Drofile 


lu  (x,y,z)  = 


.2  2 »  /  2  /  . 

a-(x  +  y  )/  a  (z) 


TT  a2(z) 


(B-5) 


The  beam  radius  a(z)  can  be  written 


a(z)  =  aQf(z)  , 


(B-6) 


where  is  the  beam  radius  in  the  source  plane  z=0.  For  the  velocity  we  assume 


v(z)  =  v  g(z) 
o 


(B-7) 


where  v  is  the  velocity  at  z=0.  Changing  to  the  normalized  variables 


u  =  x/a(z) 
w  =  y/a(z) 


(B-8) 


and  using  Eqs .  (B-5)  -  (B-7)  in  Eq.  (B-3)  gives 


=  *oP 


n  n  a" 


(x  h  +  y  J 

\  -  CO  ^ 


(B-9) 


where  e  =  e  (z=0)  is  a  constant  and 

°  ,  2  2. 

/  X  -(u  +  w  ) 
h(u,w)  =  e  v 


(B-10) 


B-2 
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Using  Eq.  (B-9)  in  Eq.  (B-2)  leads  to  the  general  expression  for  the  thermal 
distortion  exponent 


f  C  *  ('  —  +  b  -  )  +  5  (  -  +  -  -  )  3  b  . 

TT  a3  o  au  gu  h  aw  f(z) 


+  1  —  )+y(L  +  ii2 


p*  r  *  ^  -  ra  oh  e 

j  [xh+y  J  _  flu*]  _ _  dz"  iz'l 

o  *•  3w  g(z")  f^z")  I 


(B-ll) 


To  obtain  a  tractable  expression  for  the  distortion  parameter  it  is  necessary  to 
ignore_the  z-  dependence  in  h  and  the  coordinates  u^v,  so  that  the  terms  within 
the  square  brackets  can  be  removed  from  the  z-  integrations.  This  can  be  regarded 
as  being  equivalent  to  evaluating  Y  near  the  axis  (i.e.,  for  u  &  w  «  0)  where  it 
becomes  simply  a  constant  indicating  the  change  in  on-axis  intensity.  Equation 
(B-ll)  thus  becomes 


Y  =  N 
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where 
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is  the  general  expression  for  the  distortion  parameter  for  the  case  of  an  un¬ 
collimated  beam  with  a  varying  velocity  along  the  propagation  path.  N  is  the 
familiar  (e.g.,  see  Ref.  19  Eq.  (2/^/  collimated  beam  distortion  parameter  given 

by 
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which  is  valid  for  otz  «1  for  a  laser  beam  of  radius  a0  with  the  uniform  wind 
velocity  v0  .  The  term  in  brackets  in  Eq.  (B-13)  can  be  regarded  as  a  correction 
factor  that  modifies  to  account  for  finite  01  z,  focusing  and  a  variable  wind  or 
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sieving  velocity. 


The  form  of  the  correction  factor  for  focused  beam  propagation  vith  a  uniform 
vind  velocity  is  derived  first.  Thus,  taking  g(z)  =  1  and  assuming  the  linear 
beam  profile 
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vhere  a^  =  a(z^)  is  the  beam  radius  at  the  focal  range  z  ,  and  also  assuming 
orz  «  1,  the  evaluation  of  Eq.  (B-13)  is  straight  forwardfand  yields  the  result 
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vhere  x  =  a^/a^,  .  The  focused  beam  N  in  Eq.  (B-l6)  differs  from  the  incorrect 
expression  derived  previously  (Ref.  1,  Eq.  (IT)  )  on]y  by  the  factor  q  vhich  is 
plotted,  as  a  function  of  aQ/af  in  Fig.  B-l.  For  the  case  aD/af  <  1,  the  beam 
is  diverging  while  for  a0/af  >  1,  the  beam  focuses.  In  the  limit  of  very  strong 
focussing  ,  i.e.,  a  /af  ->  <*> ,  the  factor  q  approaches  2  and  hence  -»  2  N 

(a0/af)- 


Next,  the  case  of  slewing  is  considered  for  a  collimated  beam  (i.e.,  for 
f(z)  =  i)  with. .a-z  «  1.  For  this  situation  the  velocity  can  be  written 


v(z)  =  v0  +  Qz  , 

where,  as  shown  in  Fig.B-2,the  slewing  rate  is 
.1  =  (vf  -  v0)/zf  , 
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and  Vf  -  v(zf). 


From  Eqs.B-7  and  B-l8  we  have  that 
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CORRECTION  FACTOR  FOR  FOCUSED  OR  DIVERGING  BEAM  N 
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FIG.  B-2 
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Influence  of  Mirror  Aberration  on  the 
Thermal  Blooming  in  a  Wind 

by 

Frederick  G.  Gebhardt 
United  Aircraft  Research  Laboratories 
East  Hartford,  Comecticut  06108 

ABSTRACT 

Experiments  comparing  the  thermal  distortion  effects  of  wind  for  laser 
beams  with  and  without  astigmatism  are  described.  The  astigmatism,  which  is 
obtained  by  using  a  spherical  mirror  at  a  45°  angle  of  incidence,  produces 
significant  ellipticity  in  the  undistorted  beam  that,  in  addition,  varies 
continuously  along  the  propagation  path.  The  results  obtained  with  the  aberrated 
beam  show  an  ~-20$  higher  maximum  peak  intensity  as  a  function  of  laser  power 


than  the  unaberrated  beam. 
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Influence  of  Mirror  Aberration  on  the 
Thermal  Blooming  In  a  Wind 

by 

Frederick.  G.  Gebhardt 
United  Aircraft  Research  Laboratories 
Bast  Hartford,.  Connecticut  06108 

A  high  power  cw  CO2  laser  beam  propagating  in  the  atmosphere  can  be 
severely  bent,  distorted,  and  spread  by  the  self-induced  thermal  lens  effects 
of  a  transverse  wind.  Thus,  it  is  of  interest  to  examine  techniques  for 
reducing  or  minimizing  the  deleterious  effects  of  thermal  blooming  on  the 
propagation  of  CCg  laser  radiation  in  the  atmosphere.  One  approach  to  this 
problem  is  to  look  for  source  distributions  for  the  laser  transmitter  which 
experience  less  beam  degradation  than  say,  the  focused  gaussian  source.  For 
example,  numerical  calculations  have  shown  that  the  form  of  the  initial 
intensity  jrofile  has  a  significant  effect  on  the  details  of  the  thermal  distor¬ 
tion  of  the  laser  beam.  (•*■)  The  use  of  a  specially  contoured  mirror  to  modify 
the  source  wave  fiont  shape  is  also  being  considered  as  a  possible  way  of 
reducing  thermal  blooming  effects . (2 )  In  this  report  we  describe  an  experiment 
in  which  the  effect  of  a  mirror  aberration  on  the  thermal  distortion  by  a  wind 
is  investigated.  The  astigmatism  of  a  tilted  spherical  mirror  was  used  as  a 
simple  means  for  modifying  the  source  wave  front  shape.  The  objective  of  the 
experiment  was  simply  to  determine  if  the  astigmatic  phase  distortion  can 
provide  any  reduction  in  the  thermal  distortion  effects  over  the  unaberrated 
beam  case. 
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Description  of  Experiment 

The  thermal  distortion  effects  of  wind  were  examined  experimentally  by 
propagating  CO2  laser  radiation  through  liquid  CS2.  The  COg  laser  beam  was 
directed  downward  through  a  5  cm-diameter  by  12.7  cm-length  glass  cell  contair.  •• 
mg  CS2  and  equipped  with  NaCl  wondbws  at  each  end.  The  liquid  cell  was  moved' 
across  the  beam  with  a  mechanical  scanner  at  the  velocity  of  1  cm/sec  to 
simulate  the  effect  of  a  wind.  The  laser  beam  was  focused  on  the  detector, 
which  was  situated  hear  the  cell  exit  window,  with  a  spherical  copper  mirror 
of  68  cm  radius  of  curvature.  As  shown  in  Fig.  1,  the  plane  of  incidence 
(defined  by  tne  incident  laser  beam  and  the  normal  io  tne  norrory.  is  jperpendieular 
to  the  direction  of  motion  of  the  liquid  cell.  As  the  angle  of  incidence  $, 
between  the  laser  beam  and  the  normal  to  the  mirror  increases  from  zero,  the 
aberration  known  as  astigmatism  increases  in  strength. (3)  For  $  =  8°,  the 
astigmatism  is  negligibly  small  for  our  purposes,  and' we  refer  to  this  as  the 
unaberrated  case.  For  the  aberrated  case,  as  shown  in  Fig.  1,  $  =  45°,  which 
results  in  a  significant  amount  of  astigmatism. 

To  insure  that  differences  observed  in  the  thermal  distortion  between  the 
aberrated  and  unaberrated  cases  could  be  attributed  solely  to  the  astigmatism, 
special  care  was  taken  to  keep  the  distance  between  the  laser  output  mirror  and 
the  curved  focusing  mirror  the  same  for  each  case.  The  18  cm  distance  between 
the  curved  mirror  and  the  liquid  cell  was  needed  to  allow  clearance  between  the 
laser  beam  and  cell  for  the  small  (8°)  angle  of  incidence  for  the  unaberrated 


case. 


L921004-8 


The  cross  hatched  patterns  shown  in  Fig.  1  represent  the  undistorted  laser 
beam  shapes  at  the  planes  indicated.  These  were  determined  from  measured 
intensity  profiles  (figs.  3  and  4)  and  photographs  of  ihe  beam  taken  with  Kalvar 
film  (See  Ref.  1,  p.  36).  For  the  unabefrated  case,,  the  laser  beam  is  circular 
and  focuses  slightly  through  the  moving  liquid  cell.  For  the  aberrated  case, 

V 

X  / 

,  / 

on  the  other  bend,  the  astigmatism  produces  a  pronounced  ellipticity  in  the 
laser  beam  shape  at  the  cell  entrance .  The  major  axis  of  the  elliptical  pattern 
is  parallel  with  the  cell  motion  simulating  the  wind.  As  the  beam  propagates 
from  the  cell  entrance  to  the  detector  the  beam  shape  becomes  circular  and  the 
diameter  is  approximately  equal  to  the  mean  diameter  of  the  elliptical  input 

beam.  The  effect  of  astigmatism  is  to  focus  rays  lying  in  the  plane  of  incidence 

/ 

at  the  distance  (v  cos  5  )/2  (known  as  the  tangential  focus),  while  rays  lying 
in  the  plane  normal  to  thi e-  are  focused  at  the  distance  r/ (2  cos  § ) (the  saggital 
focus),  where  r  is  the  mirror  radius  of  curvature. The  observed  beam  patterns 

are  consistent  with  this  since  for  t  =  45° ,  the  minimum  beam  dimension  transverse 

/ 

to  the  wind  should  occur  at  the  tangential  focus,  which  is  approximately  midway 
through  the  cell.  Beyond  the  detector,  which  is  located  nearly  at  the  "circle 
of  least  confusion,"  the  beam  again  becomes  elliptical  but  with  the  major  axis 
normal  to  the  cell  motion.  The  maximum  ellipticity  with  this  orientation  should 
occur  about  13.5  cm  beyond  the  detector,  at  the  saggital  focus. 

Experimental  data  for  the  thermal  distortion  for  both  the  aberrated  and 
unaberrated  beams  were  obtained  by  measuring  the  intensity  profiles  through 
beam  center  along  the  "wind"  direction.  This  >?as  accomplished  by  scanning  the 
detector  across  the  beam  along  with  the  liquid  cell.  The  detector  was  equipped 
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with  a  5tyi  pin  hole  aperture  and  has  been  described'  in  Ref.  1.  The  pertinent 
experimental  parameters  are  summarized  in  Table  I  for  both  the  aberrated  and 
unab errs.ee d  beam  cases. 


Experimental  Results 

Images,  of  the  thermally  distorted  laser  beam  patterns  were  obtained  using 

•N, 

Kalvar  film  and  are  shown  in  Fig.  2  for  input  powers  varying  frem  0.5  to  7.0  W. 
The  "wind"  direct, ion  is  from  right-to-left  and  the  relative  beam  deflections 
can  be  measured  from  the  right  edges  of  the  pieces  of  film.  The  most  striking 
feature  of  the  distortion  obtained  with  the  aberrated  beam  is  the  narrow  width 
and  small  curvature  of  the  patterns  as  compared  with  the  unaberrated  case. 

The  measured  intensity  profiles  for  the  unaberrated  and  aberrated  beams 
are  shown  in  Figs.  3  and  4,  respectively.  The  "input  beam"  profile  refers  to 
the  beam  at  the  cell  entrance  and  the  remaining  profiles  were  measured  ~3  cm 
from  the  cell  exit  window.  The  various  profiles  are  labeled  with  both  the 
input  laser  beam  power,  P,  as  well  as  the  value  of  N,  which  is  evaluated  using 
the  relation: 


N  = 


(-dn/dT)2Pz 
noP  cPvai 
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-  e"0*) 

at 
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where  n0,  dn/dT,  p,  cp,  a  and  v  are,  respectively,  the  refractive  index,  the 
index  change  with  respect  to  temperature,  density,  specific  heat,  absorption 
coefficient  and  velocity  of  the  medium.  The  remaining  quantities  are  defined 
in  Table  I  where  the  material  properties  of  CSg  are  also  included.  The  factors 
H/\  and  ai/ao  are  included  to  account  for  the  ellipticity  and  focusing  of  the 
beams,  respectively.  (-0  „  _ 
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The  case  N  =  0,  in  Figs.  3  and  4,  corresponds  to  the  undistorted  profile 
obtained  without  the  absorbing  liquid  cell  in  place.  Two  scans  of  each  of  cite 
distorted  patterns  are  shown  in  Figs.  3  and  4  and  the  repeatability  of  the  data; 
is  clearly  evident.  Comparison  of  the  data  for  the  aberrated  beam  with  the 
unaberrated  case  shows  that  the  detailed  structure  in  the  distorted  patterns  is 
very  similar.  It  should  be  pointed  out  that  the  anomalously  small  profile  for 
N  =  1.2  with  the  unaberrated  beam  (Fig.  3)  is  probably  due  to  error  in  alignment 
of  the  detector  with  the  center  of  the  beam.  As  the  distortion  effects,  or 
equivalently,  the  values  of  N,  become  larger,  the  large  thermally  induced  beam 
spreading  transverse  to  the  wind  (c.f. ,  Fig.  2)  tends  to  make  this  slight 
misalignment  unimportant  as  indicated  by  the  remaining  results  in  Fig.  3. 

Data  for  the  peak  intensity  of  the  distorted  patterns  have  been  taken  from 
the  profiles  in  Figs.  3  and;  4,  and  are  shown  in  Fig.  5  as  functions  of  the  laser 
power.  The  dashed  lines  indicate  the  peak  intensity  dependence  upon  power  for 
the  linear  propagation  case,  i.e.,  in  the  absence  of  any  thermal  distortion. 

Note  that  the  linear  attenuation  has  been  included  in  these  curves  so  that  they 
are  readily  compared  with  the  values  for  the  distorted  peak  intensity.  The 
different  slopes  of  the  two  dashed  lines  simply  reflect  the  fact  that  the  spot 
at  the  detector  is  larger  for  the  aberrated  beam  than  for  the  beam  with  no 
aberration  (c.f.,  Figs.  1,  3  and  4).  With  thermal  distortion,  the  peak  intensity 
increases  with  power  until  the  beam  spreading  due  to  blooming  begins  to  increase 
in  proportion  to  the  power.  The  peaks  of  the  solid  curves  in  Fig.  5  correspond 
to  this  condition.  For  larger  power  levels  the  increase  in  thermally  induced 
beam  spreading  with  power  becomes  dominant  and  the  peak  intensity  decreases 
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monotonically.  Except  for  very  low  power  levels,  where  the  thermal  lens 
effects  are  weak,  the  aberrated  beam  peak  intensity  is  greater  than  for  the 
unaberrated  case.  In  particular,  the  maximum  peak  intensity  achieved  by  the 
aberrated  beam  is  --20$  greater  than  that  for  the  unaberrated  beam  and,  in  addi¬ 
tion,  this  maximum  occurs  at  one-half  the  power.  Thus,  according  to  Fig.  5> 
there  is  a  range  of  power  levels,  or  equivalently,  values  of  N,  for  which  the 
aberrated  beam  provides  a  greater  peak  intensity  at  the  target  than  the  beam 
with  no  aberration. 

In  Fig.  6,  the  relative  peak  intensities  (i.e.,  the  peak  intensity 
normalized  by  the  undistorted  value)  are  plotted  versus  N.  This  shows  the 
advantage  of  the  aberrated  beam  over  the  unaberrated  beam  in  terms  of  the 
amount  of  distortion  relative  to  the  respective  undistorted  beams.  The  low 
data  point  for  the  unaberrated  beam  at  N  =  1.2  is  believed  to  be  due  to  the 
detector  misalignment;  and,  as  noted  earlier,  the  alignment  error  becomes  less 
important  with  increasing  values  of  N. 

The  normalized  beam  deflection  dependence  on  N  is  shown  in  Fig;  7.  Again, 
the  aberrated  beam  is  somewhat  less  affected  by  the  thermal  distortion  than,  the 
unaberrated  beam.  The  reduced  deflection  and  spreading  of  the  aberrated  beam 
as  indicated  by  Figs.  5>  6  and  7  is  clearly  consistent  with  the  distorted 
beam  patterns  shown  in  Fig.  2.  It  is  perhaps  worthwhile  to  point  out  that  the 
use  of  the  parameter  N,  for  characterizing  the  thermal  distortion  of  the 
aberrated  beam,  is  not  particularly  helpful  if  comparison  is  to  be  made  with  a 
significantly  different  type  of  beam,  e.g.,  the  unaberrated  beam.  In  this 
instance  it  is  more  convenient  to  make  the  comparison  on  an  absolute  basis  as 
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done,  for  example,  in  Fig.  5.  The  use  of  N  in  presenting  the  thermal  distortion 
results  in  Figs.  6  and  7  indeed  shows  that  the  distortion  effects  for  the 
aberrated  and  unaberrated  beams  do  not  depend  on  N  (as  defined  in  Eg..  (1))  in 
a  universal  way. 

Although  the  aberrated  beam  offers  some  relief  from  thermal  blooming  effects 
in  these  experiments  it  is  not  clear  to  what  extent  the  elliptical  beam  shape 
associated  with  the  astigmatism  is  responsible.  Referring  to  earlier  work  with 
a  collimated  elliptical  gaussian  beam  (ignoring  diffraction  effects)  one  expects 
a  greater  amount  of  transverse  beam  spreading  if  b^/a^  <  1,  as  in  the  present 
case.  PO  This  is  consistent  with  the  nearly  straight  aberrated  beam  patterns 
(c.f.,  Fig.  2)  as  compared  with  the  more  curved  crescent  patterns  of  the  un¬ 
aberrated  beam.  With  b^/a^  <  1  it  is  also  predicted,  however,  that  the 
intensity  reduction  and  beam  deflection  is  increased  over  that  of  a  circular 
beam  of  radius  &£  with  the  same  power.  This  appears  to  be  inconsistent  with 
the  present  results.  No  doubt  diffraction  effects  (because  of  the  small 
Fresnel  numbers  involved)  and  the  complicated  combination  of  beam  ellipticity 
and  the  aberrated  wave  front  are  all  important  here  in  the  thermal  distortion 
process. 

Summary 

Experiments  comparing  the  t/iermal  distortion  effects  of  wind  for  laser 
beams  with  and  without  astigmatism  have  teen  described.  The  astigmatic  beam 
was  obtained  by  use  of  a  spherical  mirror  at  an  angle  of  incidence  of  4p°. 

The  astigmatic  mirror  aberration  produces  significant  ellipticity  in  the  un¬ 
distorted  beam  which,  in  addition,  is  continuously  varying  along  the  propagation 
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path.  Although  the  aberrated  beam  showed  an  rJZOf,  higher  maximum  peak  intensity 
as  a  function  of  laser  power  than  the  unaberrated  beam,  the  results  are  diffi¬ 
cult  to  analyze  or  explain.  Without  a  better  understanding  of  the  mechanism 
by  which  the  astigmatism  modifies  the  thermal  lens,  effects,  it  is  difficult 

to  draw  any  firm  conclusions  of  a  general  nature  from  these  experiments.  If 

/• 

further  experiments  a.long  this  line  are  to  be  carried  out,  it  is  recommended 
that  a  simpler  situation  be  considered  together  with  the  support  of  an  analytical 
effort. 

Helpful  discussions  with  Dr.  D.  C.  Smith  and  the  assistance  of 
Mr.  A.  Guardiani  with  the  experiments  are  gratefully  acknowledged. 
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TABLE  I 


EXPERIMENTAL  PARAMETERS 


12.7  cm 

l6  cm  (cell  input  to  detector) 

1  cm/sec 

2.2 

0.1  (includes  ~10^  NaCl  window  reflection  loss) 
-0.79  x  10-3  °c-l.  nQ  =  1.63; 

1.26  g/cm3;  cp  =  0.95  j/g°C 


UNABERRATED  BEAM 

ABERRATED  BEAM 

*  a^,  mm 

0.67 

0.86 

b^,  mm 

0.67 

0.38 

nun 

0.44 

0.67 

**F  =  kg  af/z 

2.5 

4.1 

n/p,  w"1 

12.2 

10.9 

Cell  length  t  = 
Total  path  length  z  = 
Velocity  v  = 

Oft  at 

Cell  transmission  T  = 
CS2  Properties:  dn/dT  = 

P  *5 


*  a^  =  input  l/e  beam  radius  -  parallel  with  wind 
b^  =  input  l/e  beam  radius  -  perpendicular;  to  wind 
a0  =  output  l/e  beam  radius 

**ke  =  (2tt/\z)  f  n(z')dz'  -  effective  propagation  constant 
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The  effect  of  kinetic  cooling  due  to  the  absorption  of  10.6-p  C02  laser  radiation  in  the  at' 
mosphere  has  recently-been  predicted.  In  this  letter,  experimental  evidence  for  this  effect 
is  presented.  Transient-density  increases  due  to  the  kinetic  cooling  have  been  observed  in 
COj-Nj  mixtures  following  the  propagation, of  an  ~  1-jisec  0. 5-J  C02  laser  pulse  through  the 
gas.  The  experiments  were  performed  with  varying  concentrations  of  C02  and  the  observed 
cooling  and  subsequent  heating  effects  are  explained  using  the  Usual  three-level  approxima¬ 
tion  for  vibrational-  energy  transfer  in  CC^-fy,  together  with  known  relaxation  rates. 


The  propagation  of  high-power  C02  laser  radiation  in 
the  atmosphere  can  lead  to  various  self-induced  thermal 
lens  effects  which  distort  or  otherwise  alter  the  beam 
itself. 1-5  These  effects  are  caused  by  the  changes  in 
temperature,  and  hence,  refractive  index  that  result 
from  the  absorption  of  10.6-p  radiation,  by  C02,  and 
water  vapor.  The  absorption  process  involves  vibra¬ 
tional  transitions  and  the  increase  in  vibrational  energy 
is  transferred  by  collisional  relaxation  processes  to 
translation  causing  the  gas  to  eventually  be  heated.  In 
the  case  of  K20,  as  with  many  other  absorbers,  the  re¬ 
laxation  processes  are  so  rapid  that  for  all  practical 
purposes  the  absorbed  energy  is  transferred  instanta¬ 
neously  into  ftet.t.  For  C02,  on  the  other  hand,  the  re¬ 
laxation  processes  are  much  slower  and  the  possibility 
of  a.  nonequilibrium  cooling  of  the  gas  by  absorption  has 
been  reported. 3,1  With  cooling,  the  gas  density  and  re¬ 
fractive  index  increase,  which  tends  to  focus  the  laser 
beam  in  contrast  with  the  usual  defocusing  or  “bloom¬ 
ing’’  effects  due  to  heating.  Although  the  cooling  is  dis¬ 
tinctly  a  transient  effect  and  is  thus  important  for  pulsed 
propagation, 4  Wallace  and  Camac5  have  shown  its  im- 
,portance  also  for  cw  propagation  with  a  wind  or  beam 
slewing. 

In  this  letter,  experimental  evidence  for  the  kinetic 
cooling  of  a  gas  by  the  absorption  of  C02  laser  radiation 
is  presented.  Transient-density  increases  have  been 
observed  in  C02-N2  mixtures  following  the  propagation 
of  an  ~  1-psec  0. 5-J  C02  laser  pulse  through  the  gas. 

The  observed  cooling  interval  and  delay  in  gas  heating 
increase  with  decreasing  C02  concentration  in  accor¬ 
dance  with  known  relaxation  rates. 

A  detailed  theoretical  treatment  of  the  energy-exchange 
processes  involved  with  10.6-p  absorption  in  air  has 


been  presented  by  Wood,  Camac,  and  Gerry. 4  For  the 
present  experiment  where  gas  mixtures  of  varying  con¬ 
centrations  of  C02  with  N2  have  been  used  rather  than 
natural  air,  the  kinetic  model  must  be  changed  some¬ 
what  from  that  in  Ref.  4,  particularly  with  regard  to 
the  rates  involved.  Also,  the  use  of  a  pulse  width 
shorter  than  the  important  relaxation  times  alters  the 
effects  in  the  experiment  somewhat  from  the  predictions 
of  Ref.  4.  In  order  to  discuss  the  absorption  of  10. 6~p 
radiation  and  subsequent  relaxation  processes  in  a  C02- 
N2  mixture,  the  familiar  three-level  approximation  for 
the  C02  laser8  in  Fig.  1  is  convenient.  Group  (0)  is  the 
ground  level  with  no  vibrational  energy;  level  (1)  in¬ 
cludes  the  bending  (i/2)  and  symmetric  stretch  (n,) 
modes  of  C02  with  the  lower  level  (100)  of  the  ’.0. 6-p 
transition;  and  the  level  (2)  includes  the  first  vibra¬ 
tional  level  of  N2  together  with  the  upper  level  (001)  of 
the  absorbing  transition,  which  is  the  asymmetric 
stretch  mode  (vz)  of  C02.  Within  the  limits  of  the  three- 
level  approximation  the  vibrational  levels  within  the  in¬ 
dividual  groups  are  assumed  to  be  closely  coupled  by 
rapid  relaxation  processes;  and,  the  only  important  re¬ 
laxation  times  are  t2!  and  t10,  which  govern  the  energy 
transfer  between  the  levels  (2)-(l)  and  (l)-(0),  respec¬ 
tively.  The  absorption  of  a  10. 6-p  C02  laser  photon 
[the  absorption  coefficient  for  7=300^  and  1-atm 
pressure  is  aaXCOz  (2X10*3)  cm*1,  where  Xc0  is  the 
mole  fraction  of  C02  (Ref.  7)]  produces  a  001  vibra- 
tionally  excited  molecule  and  eliminates  one  molecule 
in  the  100  state.  The  nearly  resonant  energy  transfer 
between  the  C02  (001)  and  N2  (i>=  1)  levels  rapidly  es¬ 
tablishes  vibrational  equilibrium  and  the  combined  en¬ 
ergy  levels  (2)  relax  through  collisions  to  (1)  at  the  rate 
t2}.  Level  (1),  which  has  been  depleted  by  the  absorp¬ 
tion  process,  is  restored  to  thermodynamic  equilibrium 
either  through  the  ground  state  (0)  at  the  rate  rjj,  or  die 


D-UL- 

Appl.  Phys.  Lett.,  Vol.  20,  No.  3,  I  February  1 972 


L9210C&-8 


F.G.  GEBHARDT  AND  D.C.  SMITH 


scope  and  used  to  illuminate  the  interferometer.  A  ger¬ 
manium  lens  L  with  25-cm  focal  length  is  used  to  focus 
the  10.6-m  radiation  into  the  gas  cell  through  an  anti¬ 
reflection-coated  germanium  window-  W  that  also  serves 
as  a  mirror  for  the  visible  beam  in  the  interferometer. 
The  density  changes  in  the  gas  cell  that  result  from  the 
absorption  of  CQ  laser  radiation  are  measured  with  a 
photodiode  placed  in  the  fringe  field  at  the  interferom¬ 
eter  output.  With  the  interferometer  adjusted  to  pro¬ 
duce  an  infinite  fringe,'  a  change  in  gas  density  modifies 
the  interferometer  intensity  and  a  signal  is  obtained 
from  the  photodiode  that  is  amplified  (Tektronix  1121, 
5-Hz  17-MHz  bandwidth)  and  displayed  on  an  oscillo¬ 
scope.  The  combined  photodiode-amplifier  rise  time 
is  <  2  Msec  as  measured  with  a  spark  gap.  Rotation  of 
the  optical  flat  F  is  used  to  bias  the  interferometer  for 
maximum  sensitivity  and  for  calibration  by  establishing 
the  magnitude  and  direction  of  the  phase  change.  The 
CQ,  laser  beam  was  measured  from  burn  patterns  on 
thermofax  paper  to  be  ~  2  x  5  mm  at  the  cell  entrance 
and  exit  windows,  and  the  focused  spot  size  midway 
through  the  cell  was  estimated  to  be  ~0. 2  mm.  The 
photodiode  aperture  was  4  mm  in  diameter. 

The  experiment  consisted  essentially  of  monitoring  the 
density  changes  in  the  CCfe-Nj  mixture  following  the  ab¬ 
sorption  of  a  portion  of  the  C02  laser  pulse  energy.  The 
experiments  were  performed  with  XCOz  varied  from 
9. 05  to  1. 0  at  a  total  pressure  of  1  atm,  and  the  results 


FIG..  1.  Simplified  vibraiional-energy-level1  diagram  showing 
the  three-level  approximation  for  COj-Nj  gas  mixtures.  - 


deactivation  of  level  (2)  at  the  rate  Tz\f  depending  on 
which  is  faster.  With  the  former  process  (i.e.,  t10<t21) 
vibrational  equilibrium  is  established  at  the  expense  of 
translational  energy,  and  the  gas  is  cooled  until  level 
(2)  can  relax.  In  the  latter  case  (i.e.,  r21<  t10)  no  cool¬ 
ing  occurs  and  the  gas  is  heated  following  the  collisional 
deactivation  of  levels  (2)  and  (1).  In  both  cases,  gas 
heating  due  to  the  absorbed  energy  occurs  on  a  time 
scale  of  ~r21+  f10.  The  relaxation  times  r2.  and  r,0  have 
been  calculated  using  the  expressions  and  rates  given 
in  Ref.  6.  The  results  are  shown  in  Fig.  2  as  functions 
of  -^co2.  the  mole  fraction  of  C02,  for  a  pressure  of  1. 
atm.  The  characteristic  heating  time  r21+  ria  is  in¬ 
cluded  in  Fig.  2  as  the  broken  curve.  Decreasing  *co2 
or  equivalently  increasing  the  N2  concentration  Xt,2  sig¬ 
nificantly  increases  the  upper  level  (2)  lifetime  r21  while 
slightly  decreasing  the  lower  level  (1)  lifetime  r10.  Ac¬ 
cording  to  Fig.  2,  the  kinetic  cooling  effect  should 
occur  for  Xcl,2  <  0. 6,  since  under  these  conditions 
r21>  t)0.  As  Xcc>2  decreases,  the  cooling  effect  becomes 
more  pronounced  since  r2I  increases.  The  magnitude  of 
the  cooling  and  subsequent  heating,  however,  depends 
on  the  absorption  coefficient  which  is  proportional  to 
XCo2.  It  should  be  pointed  out  that  the  kinetic  model 
used  for  air  in  Ref.  4  differs  somewhat  from  the  three- 
level  C02-N2  model  because  of  the  low  (~  0. 03%)  C02 
concentration  and  the  presence  cf  H20  and  02  in  air. 
Although  this  modifies  the  various  relaxation  rates  in¬ 
volved,  essentially  the  same  vibrational-energy-trans¬ 
fer  processes  in  C02  and  N2  are  responsible  for  the 
cooling  effect  in  both  cases.  The  main  difficulty  in  ob¬ 
serving  experimentally  the  cooling  effect  in  normal  air 
is,  of  course,  due  to  the  small  absorption  coefficient. 

In  addition,  while  H20  niters  '.he  relaxation  rates  in¬ 
volved  in  the  cooling  process,  it  also  absorbs  laser 
radiation  and  causes  beam-heating  effects  that  compete 
with  the  cooling. 

The  experimental  arrangement  is  shown  in  Fig.  3.  A 
1-Msec  0. 5-J  pulse  from  an  electrically  excited  C02 
laser8  is  passed  through  a  gas  cell,  12  cm  long  by  3. 8 
cm  in  diameter,  containing  a  mixture  of  C02  and  N2. 

The  gas  cell  is  equipped  with  NaCl  windows  and  is 
placed  in  one  arm  of  a  Macn-Zehnder  interferometer. 

A  He-Ne  laser  beam  at  6328  A  is  expanded  by  a  tele- 
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FIG.  2.  Relaxation  times  calculated  U3ing  data  from  Ref.  6  for 
the  three-level  approximation  hi  COj-N2  mixtures  are  shown 
together  with  experimentally  determined  thermal  rise  times. 
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KINETIC  COOLING  OF  GAS 


FIG.  3.  Experimental  arrangement  for  the  interferometric  ob¬ 
servation  of  kinetic  cooling.  L:  Ge  lens,  25-cm  focal  length; 

W;  Ge-window  Ar  coated  for  10.6  p;  B:  beamsplitter;  M:  mir¬ 
ror;  F;  optical  flat;  C:  gas  cell,  3. 8  cm  in  diam;  P:  vacuum 
gauge. 

for  JfCo2=0.1,  0.2,  0.3,  0.4,  0.6,  and  l.Oare  shown 
in  Fig.  4.  The  photodiode  signals  in  Fig.  4  show  the 
normalized  density  changes  with  time  (horizontal  scale 
=  50  psec/div),  and  a  density  increase  corresponding  to 
the  cooling  effect  is  observed  for  the  cases  where 
XCOz  <  0. 6  as  predicted.  For  Xcc>2  >  0. 6,  only  a  density 
decrease  due  to  heating  is  observed  with  essentially  no 
change  in  the  signal  as  XCot  is  increased.  The  density 
change  can  be  expressed  In  terms  of  the  photodiode- 
signal  voltage  since  it  is  small  compared  with  that  re¬ 
quired  to  change  the  optical  path  length  by  ix,  as  deter¬ 


mined  by  rotating  the  optical  flat  F.  Thus,  the  photo¬ 
diode  signal  is  assumed  to  be  nearly  proportional  to  the 
interferometer  path  difference  and,  hence,  the  gas  den¬ 
sity.  Since  the  laser  pulses  are  short  compared  with 
the  relaxation  times  shown  in  Fig.  2,  the  gas-tempera¬ 
ture  changes  following  the  laser  pulse  occur  an  the 
same  time  scale  as  these  vibrational  decay  processes. 
The  gas-density  change,  and  hence  the  interferometer 
signal,  on  the  other  hand,  are  delayed  by  the  acoustic 
transit  time  rt=D/v„  where  D  is  the  beam  diameter 
and  v,  the  sound  velocity.  Following  this  time  the  gas- 
density  change  with  temperature  for  a  constant  pressure 
is  observed.  Thus,  considering  the  case  of  pure  C02 
(i.e.,  XCOt=  1.0),  the  gas  heating  should  occur  after  r21 
+  t10ss  17  psec  (see  Fig  2),  and  the  acoustic  time  pre¬ 
ceding  the  gas  density  change  is  r4«  13  psec,  where  we 
have  assumed  D=  3. 5  mm  and  t>,  =2.8x  io4  cm/sec  for 
C02  at  21  °C.  The  ~  10—  15-psec  rise  time  of  the  inter¬ 
ferometer  signal  for  XCOi=  1.0  in  Fig.  4  is  consistent 
with  this  estimate  for  t„.  That  the  rise  time  is  asso¬ 
ciated  with  the  acoustic  time  has  been  established  by 
reducing  the  pressure  to  ~  0. 5  atm  and  observing  only 
a  small  change  in  the  signal  rise  time  in  spite  of  the 
doubling  of  the  relaxation  times  at  the  lower  pressure. 
As  the  C02  concentration  decreases  below  60%,  r10  be¬ 
comes  shorter  than  r2t,  and  the  gas  cools  after  the  time 
r10~6  psec.  The  rise  time  of  the  interferometer  signal 
showing  the  density  increase  is,  of  course,  delayed  by 
the  acoustic  transit  time.  Since  the  heating  time 
(r21  +  t,0)  increases  with  decreasing  XCOz,  becoming 
much  greater  than  the  acoustic  transit  time,  this  time 
eventually  dominates  the  time  for  the  gas  density  to  de¬ 
crease  to  the  minimum  value  after  the  cooling  interval. 
The  thermal  rise  time,  which  is  defined  as  the  time  for 


APK  =  0.4  x  10’3/Jiv  a  p\p,  ~  1.0  *  IO°/di* 
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FIG.  4,  Photodiode  signals  showing  inler- 
fcrometrically  determined  density  changes 
in  COj-Nj  gas  mixtures  following  the  prop¬ 
agation  of  an  ~l-psec  0. 5-J  C02  laser 
pulse. 
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the  density  to  reach  90%  of  the  total  change  from  maxi¬ 
mum  to  minimum  density,  has  been  measured  from  the 
oscilloscope  traces  and  the  results  are  shown  in  fig.  2. 
The  beginning  of  the  time  interval  is  taken  as  the  break 
in  the  oscilloscope  trace  which  is  caused  by  an  electri¬ 
cal  noise  transient  associated  with  the  spark  gap  used 
to  trigger  the  COz  laser.  The  measured  thermal  rise  or 
heating  time  is  seen  to  be  in  close  agreement  with  the 
calculated  heating  time  (t2,+  t19).  it  should  be  pointed 
out,  however,  that  this  agreement  is  somewhat  fortu¬ 
itous  and  probably  due  to  the  effect  of  the  finite  acoustic 
transit  time.  For  the  case  where  the  acoustic  time 
and  r10  become  negi  gible  compared  with  tZI)  the  so- 
called  thermal  rise  time  should  in  theory  approach 
~2r21.  Also,  in  this  limit,  the  density  increase  for  very 
short  times  (t«  r21)  should  be  1.44  times  the  density 
decrease  for  / »  r2I,  which  is  the  ratio  of  the  100  vibra¬ 
tional  energy  to  the  photon  energy  hv  10  #. 

For  very  long  times  the  decrease  in  gas  density  must 
eventually  decay  to  zero,  due  to  thermal  conduction. 

The  thermal-conduction  time  rc  of  ~  50  msec  has  been 
measured  for  pure  C02,  and  since  tc=D3/16x,  where 
X=  0. 11  cm2/sec  Is  the  thermal  diffusivity  for  C02,  we 
obtain  D  «  3  mm  for  the  effective  beam  diameter  within 
the  gas  cell. 

It  is  interesting  to  note  that  the  density  change  with 
A'C02=  0. 1  is  roughly  about  one-half  that  obtained  with 
pure  C02  rather  than  one-tenth  as  expected.  This  indi¬ 
cates  that  the  C02  laser  pulse  energy  density  exceeds 
the  saturation  value  for  pure  C02,  and  when  XCo2  is  re¬ 
duced  by  adding  N2  the  saturation  energy  density  in¬ 
creases  to  compensate  somewhat  for  the  reduction  in 
thejabsorption  coefficient.  For  a  laser  pulse  shorter 
than  the  upper-level  relaxation  processes,  as  in  the 
present  case,  the  absorbed  energy  density  under 
strongly  saturated  conditions  is  a0Es  =  nlhui0tt/2,  where 
«i  is  the  number  density  of  G02  molecules  in  the  lower 
absorbing  (100)  level  and  <*„  is  the  unsaturated  absorp¬ 
tion  coefficient.  That  is  to  say,  E,  is  the  energy  per 
unit  area  required  for  one-half  of  the  (100)  C02  mole¬ 
cules  to  be  excited  by  absorption  to  the  (001)  leveLand 
thus  approximately  saturate  the  transition.  For  pure 
C02  at  T=  300  °K  and  1-atm  pressure. with  a0=  2x  10° 
cm'1  and  n,=  7. 25*  1016  cm'3,  we  obtain  for  E,  the  value 
0. 34  J/cm2.  The  laser  pulse  enercy  of  0. 5  J  together 
with  thebeam  radius  of  1. 5  mm  (based  on  the  measured 
thermal-conductioh  signal  decay)  gives  ~7  J/cm2,  which 
clearly  exceeds  the  estimated  saturation  energy  density. 


The  observed  density  decrease  in  Fig.  4  for  pure  C02 
corresponds  to  a  temperature  rise  of  AT~0. 5cC.  As¬ 
suming  saturated  conditions  and  using  the  energy  equa¬ 
tion,  we  find  that  £*-0. 38  J/cm2  for  pure  CCh,  which 
is  in  good  agreement  with  the -calculated  value. 

It  should  also  be  noted  that  the  small  pulses  appearing 
in  the  interferometer  signals  in  Fig.  4  are  associated 
with  repeated  reflections  of  the  acoustic  wave  from  the 
walls  of  the  gas  cell.  Their  occurrence  in  pure  C02  at 
intervals  of  ~  125  psec  after  the  initial  density  change 
corresponds  well  with  the  3.8-cm  tube  diameter.  For 
the  case  of  XCOj=  the  pulse  interval  is  reduced  to 
about  100  psec,  which  can  be  explained  by  the  increase 
in  sound  velocity  to  vs~  3. 4x  104  cm/sec  due  to  the  ad¬ 
dition  of  N2. 

In  summary,  these  results  have  shown  the  first  experiv 
mental  evidence  of  kinetic  cooling  by  absorption  of  laser 
radiation.  The  cooling  and  subsequent  heating  occur  on 
the  time  scale  of  known  vibrational  relaxation  rates, 
and  the  simple  three-level  model  for  the  C02-N2  system 
appears  to  be  adequate  for  modeling  the  effect.  In  the1 
■case  of  atmospheric  propagation,  the  times  associated 
with  the  cooling  will  be  increased  to  milliseconds  or 
greater  because  of  the  longer  lifetime,  t2J  associated 
with  the  lower  C  (^concentration.  The  kinetic  cooling 
can  be  a  beneficial  effect  for  propagating  pulsed  beams 
which  have  a  pulse  width  comparable  to  r21,  or  for  cw 
beams  propagating  with  a  cross  wind  such  that  the  tran¬ 
sit  time  across  the  beam  is  comparable  to  r21. 
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-  .APPENDIX  E 

Comparison  of  Numerical  and  Experimental 
Thermal  Distortion  Results 

In  this  appendix  a  measured  thermally  distorted  laser  beam  pattern  is  compared 
with  results  obtained  numerically.  In  Fig.  E-l,  the  distorted  CO5  beam  pattern 
obtained  with  Kalvar  film  in  a  liquid  CS2  experiment  simulating  a  wind  is  shown  on 
the  right-hand  side.  The  value  of  N  as  computed  using  Eq.  (1)  and  the  appropriate 
experimental  conditions  (Ref.  1,  Figs.  36  and  37,  P=!2  W  case)  is  20  for  the  measured 
pattern.  On  the  left-hand  side  of  Fig.  E-l,  numerical  results  obtained  by  Bradley 
and  Herrmann  (Ref.  2)  with  the  Lincoln  Laboratories  nonlinear  propagation  code  are 
shown  for  a  case  with  N  =  16,  which  is  roughly  comparable  to  the  distortion  parameter 
for  the  liquid  experiment'.  In  each  case  the  beams  were  focused  and  a  close  similarity 
between  the  numerical  and  experimental  results  can  be  seen,  particularly  with  regard 
to  the  overall  shape  of  the  patterns  and  the  presence  of  the  small  central  lobes. 

In  particular,  both  the  values  of  IR  and  Ax/af,  which  are,  respectively,  the  peak 
intensity  normalized  by  the  undistorted  intensity  and  the  normalized  beam  shift  into 
the  wind,  are  in  agreement  as  shown  in  Fig.  E-l.  The  calculated  constant  intensity 
contours  are  labeled  to  indicate  the  values  of  intensity  normalized  by  the  peak. 

The  constant  intensity  contour  represented  by  the  edge  of  the  measured  pattern 
corresponds  to  14  percent  of  the  peak  which,  is  indicated  by  the  x .  This  was 
determined  by  comparing  the  measured  intensity  profile  (Ref.  1,  Fig.  37,  P  =  2  W 
case)  with  the  Kalvar  image  of  the  beam  pattern.  The  dashed  circles  represent  the 
undistorted  l/e  intensity  contours  of  the  focused  spots  for  each  case. 

The  thermal  distortion  results  are  shown  in  Fig.  E-l  in  terms  of  normalized 
dimensions  and  intensity  levels;  and  the  physical  parameters  for  each  case  are  also 
given  in  normalized  fashion  in  terms  of  the  distortion  parameter  N,  the  fractional 
absorption,  az,  and  the  ratio  of  initial  to  focused  spot  sizes,  a^/a^  of  the  beams 
in  the  absence  of  distortion.  Of  these  parameters,  the  values  for  az  differ  the 
most  between  the  numerical  and  experimental  cases.  In  view  of  the  results  showing 
the  az  dependence  of  the  distortion  in  Sec.  2.4,  the  larger  az  in  the  liquid  experi¬ 
ment  may  be  responsible  for  the  good  agreement  between  the  numerical  and  experimental 
results  in  spite  of  the  fact  that  the  values  of  N  differ  between  the  numerical 
(N  =  16)  and  experimental  (N  =  20)  cases. 

Perhaps  the  most  important  point  to  be  made  by  the  comparison  in  Fig.  E-l  of 
the  numerical  results  with  those  obtained  experimentally  concerns  the  vastly  different 
conditions  between  (a)  the  laboratory  simulation  of  wind  with  a  small  (~  5  in.  long) 
cell  of  liquid  CSg  moving  across  the  beam,  and  (b)  the  numerical  calculations  which 
were  performed  using  conditions  typical  for  high-power  propagation  in  the  atmosphere. 
The  actual  conditions  associated  with  the  liquid  CS2  experiment  and  also  for  which  the 
numerical  result  in  Fig.  E-l  was  obtained  are  given  in  Table  E-I.  Thus,  from  the 
rather  good  agreement  indicated  by  the  results  in  Fig.  E-l,  there  should  be  no 
question  about  the  validity  of  the  liquid  simulation  of  atmospheric  propagation  in 
a  wind.  In  addition,  Fig.  E-l  provides  convincing  evidence  for  the  usefulness  of 
the  parameter,  N,  for  the  scaling  or  comparison  of  the  convective  thermal  lens 
effects  under  widely  differing  conditions,  e.g.,  as  in  Table  E-I. 
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TABLE  E-I 


COMPARISON  OF  CONDITIONS1  FOR  FIG.  E-l 
X  =  10.6  p, 


Numerical  Results 
(Lincoln  Labs . ) 

Liquid  CS2* 
Experiments 

Power,  P 

100  kW 

2  W 

Source  radius,  a. 

t 

10  cm 

1  mm 

Focused  radius,  af 

3.6  cm 

0.22  mm 

Range,  z 

2  km 

15.5  cm 

Wind  Velocity,  V 

2  m/s 

1  cm/s 

N 

16 

20 

0.134 

2.2 

kn  a'r/z 

0  b 

3 

5.7 

a,/af 

2.;8 

4.5  ' 

*  Liquid  CS2  results  have  been  taken  from  Ref.  1,  Figs.  36  and  37*  To 
avoid  confusion  the  distinctions  between  the  cell  length,  t  and  z, 
and  the  effective  propagation  constant,  ke  and  knQ,  are  omitted  here. 
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